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t— i Abstract: A key question in the study of M = 2 supersymmetric string or field theories is 

i — | to understand the decay of BPS bound states across walls of marginal stability in the space 
^ of parameters or vacua. By representing the potentially unstable bound states as multi- 

centered black hole solutions in M = 2 supergravity, we provide two fully general and explicit 
formula? for the change in the (refined) index across the wall. The first, "Higgs branch" 
formula relies on Reineke's results for invariants of quivers without oriented loops, specialized 
CO to the Abelian case. The second, "Coulomb branch" formula results from evaluating the 

qq symplectic volume of the classical phase space of multi-centered solutions by localization. 

We provide extensive evidence that these new formulae agree with each other and with the 
h mathematical results of Kontsevich and Soibelman (KS) and Joyce and Song (JS). The main 

t— i physical insight behind our results is that the Bose-Fermi statistics of individual black holes 

participating in the bound state can be traded for Maxwell-Boltzmann statistics, provided 
t—i the (integer) index ^(7) of the internal degrees of freedom carried by each black hole is 

replaced by an effective (rational) index £2(7) = Yl m h Q{l/ m )/ m2 - A similar map also exists 
^ for the refined index. This observation provides a physical rationale for the appearance of 

the rational Donaldson-Thomas invariant ^(7) in the works of KS and JS. The simplicity of 
the wall-crossing formula for rational invariants allows us to generalize the "semi-primitive 
wall-crossing formula" to arbitrary decays of the type 7 —> M71 + with M = 2,3. 
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1. Introduction and summary 

In quantum field theories and string theory vacua with extended supersymmetry, the spec- 
trum of BPS states can sometimes be determined exactly in a weakly coupled region of the 
space of parameters (or vacua). In extrapolating the BPS spectrum to another point in 
parameter space, one must be wary of two issues: BPS states may pair up and disappear, 
and single particle states may decay into the continuum of multi-particle states. The first 
issue can be evaded by considering a suitable index fi(7; t a ), where 7 is the vector of electric 
and magnetic charges carried by the state and t a parametrizes the value of the couplings (or 
moduli), designed such that contributions from long multiplets cancel. The index ^(7; t a ) is 
then a piecewise constant function of the parameters t a . To deal with the second problem, it 
is important to understand how ^(7; t a ) changes across certain codimension-one subspaces of 
the parameter space, known as walls of marginal stability, where a single-particle BPS state 
becomes marginally unstable against decay into two (or more) BPS states [1, 2, 3, 4, 5]. 

Initial progress in this direction for four-dimensional string vacua came from supergrav- 
ity, where BPS states are represented by (in general multi-centered) classical black hole 
solutions. Since the class of multi-centered solutions potentially unstable at a certain wall 
of marginal stability exists only on one side of the wall [6, 7, 8], the discontinuity Af2(7) in 
Jl(7, t a ) is equal to the index of the multi-centered solutions with total charge 7, up to a sign 
depending whether one enters or leaves the side on which these solutions exist [9] . Based on 
this physical picture, one easily finds that the jump of the index in the simplest case, where 
the only configuration that may appear or disappear across the wall is a two-centered solu- 
tion with primitive charge vectors 71, 72, is given by the "primitive wall-crossing formula" 
[9]: 

Afi( 7 ) = (-1)^ +1 | 7l2 | fi( 7l ) ft( 72 ) . (1.1) 

With some more effort one can also compute Afl( 7 ) in "semi-primitive" cases, where the 
relevant multi-centered solutions which appear or disappear across the wall of marginal 
stability are halos of black holes with charges in multiple of 72, orbiting around a core of 
charge 71 [9]. 

While applying this method to the general "non-primitive" case seemed out of reach up 
until now, the breakthrough came from the mathematical front, with the works of Kontsevich 
and Soibelman (KS) [10, 11] and Joyce and Song (JS) [12, 13, 14]. In these works, general 
formulae were derived for the discontinuity of generalized Donaldson-Thomas (DT) invariants 
under changes of stability conditions. It is generally believed that generalized DT invariants 
are the appropriate mathematical embodiment of the physical BPS invariants. Although the 
KS and JS wall-crossing formula? look very different, there is by now much evidence that they 
are equivalent 1 . Since these two formula? appeared, much efforts have been devoted towards 
interpreting, deriving and checking these wall-crossing formula? in various physical settings 
[15, 16, 17, 18, 19, 20, 21, 22, 23]. Our goal in this paper is to rederive the wall-crossing 
formula using multi-centered black hole solutions in supergravity and extend the original 

■""We have been informed by D. Joyce of a general argument showing the equivalence of the KS and JS 
wall-crossing formulae. The equality of the DT invariants as defined by KS and JS seems, on the other hand, 
less firmly established. 
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black hole halo picture to the general "non-primitive" case. We also carry out extensive 
comparisons between our formula? and those of KS and JS. 

One intriguing aspect of the KS and JS wall-crossing formulae is the appearance of two 
types of BPS-invariants, an integer-valued invariant Q (roughly Tr (— 1) F ) and a rational- 
valued invariant Q. The two are related by the "multi-cover formula" 



where the sum runs over all positive integers m such that 7/m lies in the charge lattice. We 
shall take (1.2) as the definition of Cl. Similar divisor sums appear in various instances in 
quantum field theory, e.g. in Schwinger's computation of pair creation in an external electric 
field (see e.g. Eqs. (4.118) and (4.119) in [24]), and in enumerative geometry, e.g. the multi- 
cover formulae for Gromov-Witten invariants [25], which indeed are naturally understood 
from pair creation arguments [26, 27]. The rational invariants Cl also arose in constructions 
of modular invariant black hole partition functions consistent with wall-crossing [28, 29, 30], 
and in studies of D-brane instanton corrections to the hypermultiplet moduli space metric 



For the purposes of computing the jump in the BPS spectrum, the rational invariants 
turn out to be especially convenient. Indeed, one consequence of the JS/KS wall-crossing 
formulae is that the variation Af2(7) of the rational invariants across a wall of marginal 
stability, when expressed in terms of the rational invariants ^(7') on one side of the wall, 
involves only "charge conserving" terms, i.e. sums of products of fi(7i) for different 7, such 
that Ylili = 7- I n contrast, the variation of the integer invariants AQi^y), expressed in 
terms of the integer invariants ^(7), does not satisfy this property, and as a result, involves 
considerably more terms. Needless to say, physical charge is conserved no matter what 
invariant one chooses to consider. 

Our main new insight, at the basis of the results presented below, is the following physical 
explanation of this phenomenon. In computing AQfr) from the index associated with a n- 
centered black hole solution carrying total charge 7 = ^™ =1 7i, each center must be treated 
as a point-like particle carrying |fi(7«)l internal states. When some of the 7, coincide, the 
corresponding centers must in addition obey Bose or Fermi statistics, depending on the sign 
of f2(7i). As a result of (anti-)symmetrizing the many-body wave-function, the total index 
associated with such a configuration involves, in addition to the product Y17=i ^(7*)> t erms 
of lower degree in f2(7«) - e.g. two identical bosons of degeneracy Q will give a degeneracy of 
Q(Q + l)/2. The terms of lower degree in Q violate the charge conservation property defined 
above. However, due to special properties of the interactions between centers (namely, the 
no- force condition between centers with mutually local charges), we show that it is possible to 
map the problem of computing the index of multi-centered black holes with individual centers 
satisfying Bose/Fermi statistics to an equivalent problem where the centers satisfy instead 
Maxwell-Boltzmann statistics. In this Boltzmannian reformulation, each center carries an 
effective (in general non integer) index ^(7) related to ^(7) via (1.2), and charge conservation 
is manifest. This provides a physical rationale of the charge conservation property of the 
wall-crossing formula written in terms of the rational invariants 0. 




(1.2) 



m|7 



[17, 31, 32]. 
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The same argument generalizes for the refined 'index' 0^(7, y), defined roughly as 
Tr (— 1) F y 2J ' A , which keeps track of the angular momentum of the BPS states. However this 
refined 'index' is only protected (i.e. immune to contributions of long multiplets) in the 
presence of a SU(2)r symmetry [21]. Such a symmetry exists in M = 2 supersymmetric 
field theory, but not in string theory or supergravity. As a result, this refined 'index' will 
in general be different at weak and strong coupling (more generally its value will depend 
on both vector multiplets and hypermultiplets). Nevertheless, one may still investigate the 
variation of Q re f(7, y) across lines of marginal stability in vector multiplet moduli space. In 
fact, KS have provided a wall-crossing formula for motivic Donaldson-Thomas invariants, 
which are conjectured to be equal to the refined invariants fi re f(7, y) at weak coupling, 
where the derived category description of D-branes is appropriate [33, 34]. Similarly, one 
may ask about the wall-crossing formula in the strong coupling region where the supergravity 
picture is appropriate. As for the standard index, we find that the variation Afi* f (7, y) can 
be computed by treating centers as Boltzmannian particles carrying internal states with 
effective refined index 



In this formulation, charge conservation is again manifest. At y = 1, f2 re f(7,?/) reduces to 



While the arguments above rely on representing BPS states as multi-centered solutions 
in supergravity, it is clear that it extends to the case of M = 2 supersymmetric gauge 
theories which can be obtained as rigid limits of supergravity theories [35]. In general, we 
expect that BPS solitons in the Coulomb phase can be represented as classical multi-centered 
solutions of the Abelian gauge theory at low energy, albeit singular ones. For the purposes 
of computing the wall-crossing, the singularity is irrelevant, and the problem can still be 
reduced to the quantum mechanics of point-like particles interacting by Coulomb law and 
scalar interactions. In particular, solitons with mutually local charges ((71,72) = 0) do not 
interact, and the above Bose-Fermi/Boltzmann equivalence carries over. 

Finally, it is worth pointing out that a similar phenomenon occurs for non-primitive 
wall-crossing in M = 4 supersymmetric string theories [36]. In this case only two-centered 
configurations contribute, and the only non-trivial effect comes from symmetrization [37, 38]. 
The variation of the index is thus given by the primitive wall-crossing formula (1.1), provided 
0(7) is replaced by with ^(7) = Ylm^^l l m ) i R this formula. Note that in contrast to 
the effective index (1.2) relevant for Af = 2 BPS states, the effective index ^(7) relevant for 
M = 4 dyons does not include any factor of 1/m 2 in its definition. This difference can be 
traced to the presence of extra fermion zero modes carried by a quarter BPS dyon in M = 4 
supersymmetric theories. Trace over these fermion zero modes for a system of identical 
particles produces an extra factor of m 2 in (1.2) compared to that for half BPS dyons in 
M = 2 supersymmetric theories. 

We shall now summarise our main results. Consider a wall of marginal stability on which 
the central charges Z 7l and Z 12 of two charge vectors 71 and 72 align. Assume further that, 
possibly after a change of basis in the lattice spanned by 71 and 72, BPS states carrying 
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fi(7) and (1.3) to (1.2). 
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charge M71 + N^ 2 exist only for (M > 0, N > 0) and (M < 0, N < 0). Then on one side 
of the wall, which we call the chamber c~, we have (71, 72) Im(Z 7l Z 72 ) > 0, and there are 
multi-centered bound states with individual centers carrying charges of the form rriiji + n i r y 2 
with different integers m^nj > 0. Here (71,72) is the symplectic inner product between 71 
and 72. On the other side of the wall, called the chamber c + , there are no bound states of this 
form. Let us denote by Q ± (a) the index Tr '(— l) 2Js on the two sides of the wall for a charge 
vector a = M71 + N~f 2 with M,N > 0. (Tr' denotes the trace after removing the fermion 
zero modes associated with broken supersymmetries.) Then the physical reasoning outlined 
above shows that the wall-crossing formula, expressed in terms of the rational invariants 
(1.2), must take the form 

n>2 {ai,...,a n } 
7=oiH ha„ 

where the sum runs over all possible unordered decompositions of a into vectors a%, . . . , a n , 
each of which is a linear combination of 71 and 72 with non-negative integer coefficients. Here, 
|Aut({«i})| is the symmetry factor appropriate for Maxwell-Boltzmann statistics, namely the 
order of the subgroup of the permutation group of n elements which preserves the ordered set 
(«i, . . . , a n ), for a fixed (arbitrary) choice of ordering. 2 Of course, one could instead decide to 
absorb this symmetry factor in the normalization of g{{a{\). The point of the normalization 
chosen in (1.4) is that g(a±, . . . , a n ) can now be identified as the index associated with 
an n-centered black hole configuration in supergravity, with the individual centers carrying 
charges a\, a 2 , . . . ,a n and treated as distinguishable particles, and furthermore carrying no 
internal degeneracy. Clearly, the same considerations imply that the wall-crossing formula 
for refined invariants takes an analog form 

n>2 {ai,...,a n } 
7=«iH ho„ 

where g re {({a>i}, y) computes the refined 'index' of the same n-centered black hole configu- 
ration, and reduces to g({cti}) at y = 1. In order to complete the wall-crossing formula we 
need to specify the factor (7ref({«j}? y) ( or its y = 1 limit g({ai}) ). While these factors can 
be extracted from the KS and JS formulae, we shall present two novel ways for computing 
them, which we call the "Higgs branch" and the "Coulomb branch" formulae. We have 
checked in many cases the equivalence of these prescriptions with the KS and JS formulae, 
although we have not yet been able to prove the equivalence rigorously. 

The "Higgs branch" formula is based on Denef's observation [7] that the spectrum of 
supersymmetric bound states of multi-centered black holes can be computed in the frame- 
work of quiver quantum mechanics. This description is appropriate at weak coupling, the 
arrows of the quiver describing the open strings stretched between two D-branes. Due to 
the fact that the charges carried by the various centers lie on a two-dimensional sublattice 
of the full charge lattice, the relevant quiver turns out to have no oriented loops. A formula 

2 Thus if the set {c^} consists of mi copies of j3i, m 2 copies of /?2 etc. then |Aut({aj})| = \\ k m^-. 
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for the motivic invariants of such quivers was given by Reineke in [39]. Furthermore, since 
the constituents of the bound states are to be treated as distinguishable particles without 
internal multiplicity, the relevant quiver carries dimension one vector spaces at each node 
(equivalently, corresponds to a U(l) n gauge theory). Reineke's formula simplifies in this 
case, leading to 

g Te f( ai ,. . .,a n ,y) = (-y)^-^^ (y 2 - l) 1 ^ £ (_i)-y E^E^^-f^f . 

partitions 

(1.6) 

Here we have denoted by = {a ir <x,) the symplectic inner product between the vectors a« 
and Oj, and have ordered the a^s such that (ai,aj) > for i < j (assuming that none of 
the vectors aj coincide). The sum runs over all ordered partitions of («i + • • ■ + a n ) into s 
vectors [3^ (1 < a < s, 1 < s < n) such that 

1. £ a 0(») = ai + ••• + «„ 

2. (3^ = Yli m i^ a i with rrif* = or 1 for each a, i. 

3. (j2 b a=l P {a \ «1 + • • • + ««) > V b With 1 < 6 < 3 - 1 

When some of the a^s coincide, the value of g Te f{cei, . . . ,a n ,y) can still be obtained from 
(1.6) by taking the limit ctij — > af even though the a^'s are supposed to be valued in the 
two dimensional lattice spanned by 71 and 72, eq. (3.1) defines a continuous function of the 
ccj's and this limit is well-defined. We have checked agreement with the KS and JS formulae 
for distinct a\ with n < 5, and in many cases where some of the a\ coincide. While it is 
not surprising that the Reineke formula is consistent with the JS formula (since the latter 
applies to moduli space of quiver representations), it is remarkable that Abelian quivers (i.e. 
quivers carrying a dimension-one vector space at each node) encode the complete information 
about wall-crossing. In Appendix §D we show that the index of certain non- Abelian quivers 
without oriented loops can be reduced to the Abelian case using the same black hole halo 
picture. 

Our second way of computing g rc f ({««}, y) instead follows from quantizing the moduli 
space of multi-centered BPS solutions, as first proposed in [40]. This description is most 
appropriate at strong coupling, when the classical supergravity description is valid. Using the 
fact that the moduli space M. n of n-centered solutions carries a natural symplectic structure 
and a Hamiltonian action of 577(2), we evaluate the integral of y 2j3 over this classical phase 
space by localization. The fixed points of the action of J3 on M. n are collinear configurations 
where all centers lie on the z-axis, with relative distances determined by 

n n 

E - a<T{ l a[J) - = A E > (i-7) 

where A is a positive real constant which can be removed by rescaling the Zj's, and a is the 
permutation which determines the order of the centers along the axis, z a u\ < z a u\ if z < j. 
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In this way we arrive at the "Coulomb branch formula" 

g mi ({«,}, y) = ^+ n -\y - y-y- n J2 s(a) y^<> , (1.8) 

permutations a 

where the sum runs over the subset of the permutation group in n elements for which the 
equations (1.7) admit a solution. The factor s(er) originates from the determinant of the 
action of J3 on the tangent space at the fixed points, and evaluates to the sign 

s(a) = (_l)#{v*(<+i)<"0} . (1.9) 

While eq. (1.8) is still implicit, since it requires solving the equations (1.7) (or rather, 
characterizing the permutations a such that (1.7) admits a solution), it provides us with an 
economic way of determining g Te { ({cti}, y), since each permutation has a unique y dependence 
and hence there are no cancellations between different permutations. In contrast both the JS 
formula and "Higgs branch" formula (1.6) involves extensive cancellations between different 
terms. We shall in fact see in §3.3 that motivated by the Coulomb branch formula one 
can find an algorithm to identify the uncancelled contributions in the Higgs branch formula 
without having to evaluate all the terms given in (1.6). 

The use of the rational invariants Cl also allows us to use the KS formula to derive explicit 
formulas for the change in the index in some special cases. This includes sectors carrying 
charges of the form 271 + N72 and 371 + for primitive vectors 71, 72 and arbitrary 
integer N. This generalizes the semi-primitive wall-crossing formula of [9] which describes 
the change in the index in the sector with charge 71 + N r j2, and some earlier results on 
higher-rank ADHM and DT invariants [41, 42, 43]. 

The rest of the paper is organised as follows. In §2 we describe how the problem of 
computing the index of multi-centered black holes can be mapped to an equivalent problem 
with the individual centers obeying Maxwell-Boltzmann statistics. We use this to derive some 
general properties of the wall-crossing formula e.g. charge conservation, and also reproduce 
the primitive and semi-primitive wall crossing formula. We also generalize the results to 
the case of refined index. In §3 we compute the index associated to n-centered black hole 
configurations in two different ways, first by mapping the problem to a quiver quantum 
mechanics and second by quantizing the classical phase space of multi-centered solutions. 
This leads to the Coulomb and Higgs branch formulae described above. In §4 we review the 
KS wall-crossing formula, and recast it in terms of the rational invariants 0, making the 
charge conservation property manifest. We verify that the motivic KS formula agrees with 
the results of §2 and §3 in many cases, and obtain higher order generalizations of the semi- 
primitive wall-crossing formula. In §5 we review the wall-crossing formula due to Joyce and 
Song and compare it to the KS, Higgs branch and Coulomb formulae. We find agreement 
in all cases that we consider. In Appendix A we illustrate the general wall crossing formula? 
in some special cases. In Appendix B we apply the results of §4 to analyze the spectrum of 
D0-D6 bound states on a Calabi-Yau 3-fold. In Appendix C we check our results against 
the spectrum of BPS states in Af = 2 SU(2) gauge theories. In Appendix D we show 
how our Boltzmann gas picture allows one to express the Poincare polynomial associated 
with quiver quantum mechanics with U(N) factors in terms of the Poincare polynomial of 
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Abelian quivers. This lends further support to the validity of the Boltzmann gas picture of 
multi-centered black holes. 

2. Boltzmannian view of the wall-crossing 

2.1 BPS states in J\f = 2 supergravity 

We consider Af = 2 supergravity in 4 dimensions, coupled to n v vector multiplets. Let M. be 
the vector multiplet moduli space, parametrized by complex scalar fields t a , a = 1, . . . ,n v , 
and T be the lattice of electromagnetic charges. T is a lattice of rank dimension 2n v + 2, 
equipped with an integer symplectic pairing. We choose a Lagrangian decomposition T = 
T m © T e , such that the symplectic pairing is given by 

(7,7') = 9AP ,A -^AeZ (2.1) 

where 7 = Go a ,<?a), l' = (p' A ,q'a)- The mass of a BPS state with charge 7 is equal to the 
absolute value of the central charge Z 7 , defined by 

Z, = e K / 2 (q A X A -p A F A ) , (2.2) 

where /C, X A and F A are appropriate functions of the moduli fields {t a }. Let "H 7 (t a ) be the 
Hilbert space of states with charge 7 in the vacuum where the scalars asymptote to t a at 
spatial infinity. The index 

n(7;*°) = -\^n^){-i) 2J *{zh) 2 (2-3) 

defines an integer £7(7; t a ), which counts the number of BPS states with sign (the overall 
coefficient ensures that a half-hypermultiplet contributes one unit to £1(7; t a )). Alternatively 
we could define the index as 

n(r,t a )= Tr^ (to) (-l) 2 ^ (2.4) 

where TV' denotes trace over BPS states, computed after removing the contribution from 
the fermion zero modes associated with the supersymmetries broken by the state. Mathe- 
matically, the BPS invariants {£7(7; t a ), 7 e T} are the generalized Donaldson-Thomas (DT) 
invariants for the derived category of coherent sheaves (in type IIA on a Calabi-Yau threefold 
X) or the Fukaya category (in type IIB on a Calabi-Yau threefold X). 

2.2 Wall crossing: Preliminaries 

The BPS invariants £^(7; t a ) are locally constant functions of t a , but may jump on codimen- 
sion one subspaces of M. (line of marginal stability), where 7 can be written as the sum 
7 = M71 + ^72, where M, N are two non-zero integers and 71 and 72 are two primitive 
(non-zero) vectors such that the phases of Z 7l and Z 72 are aligned. We denote the hyper- 
plane where the phases of 71,72 are aligned by ^(71,72)- Since the dependence of Q( / -f;t a ) 
on {t a } is mild, we shall drop t a from the argument of Q and use superscripts ± to indicate 
which side of P(7i,72) we are computing the index. 
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Clearly, ^(71,72) depends only on the two-plane spanned by 71 and 72 inside T. For 
a given choice of 7 and of this two plane, (M,N) and (71,72) are uniquely defined up to a 
common action of SL(2, Z). We shall now make a special choice of (71,72) such that there 
are no BPS states carrying charges of the form Mji — N r y 2 or —M^\ + ./V72 with M, N > 
[22]. For this reason it will be convenient to introduce the notation 



and focus on BPS states carrying charge in T. For definiteness we shall choose (71,72) < 0. 

We shall be considering the jump in the index Vt(M^i + N^) across the wall "^(71,72) 
for non-negative integers M, N. We shall denote by Z 7 the central charge associated with the 
charge 7, and choose Q + to describe the index in the chamber in which arg(Z 71 ) > arg(Z 72 ). 
In this case arg(Z A ./ 7l+A r 72 ) is greater (less) than arg(Z A: p 7l+A r/ 72 ) if M/N is greater (less) 
than M'/N'. We shall denote this chamber by c+. For Q~ the ordering of arg(Zjtf 7l+ jv 72 ) is 
reversed, and the corresponding chamber will be called c_. Using the fact that (71,72) < 
we now get in the chamber c + 



This is the condition under which two-centered bound states of black holes carrying charges 
M71 + ./V72 and M'71 + iV'72 do not exist [6]. Thus f2 + 's label the index in the cham- 
ber in which there are no two centered black holes with each center carrying charge in 
T. Generalizing this argument (see §3.2) one can show that in the chamber c + there are 
also no multi-centered black holes carrying charges of the form (M^i + N^) for different 
({M{}, {Ni}). In contrast f2~'s label the index in the chamber where such bound states are 
present. 

Note that fi + can receive contribution both from single and multi-centered black holes, 
but these multi-centered configurations consist of centers whose charges lie outside the sub- 
lattice spanned by 71 and 72. Due to this the distances between the centers remain finite 
as the wall ^(71,72) is approached. In contrast the distance between the centers carrying 
charges in T - appearing in a configuration contributing to Q~ - goes to infinity in this limit. 
Thus the configurations which contribute to Q + , even though not necessarily single centered 
black holes, can be treated as a single unit near this wall of marginal stability. For this 
reason we shall refer to Q + (Mji + Nj 2 ) as the index associated with a black hole molecule of 
charge M71 + ./V72. Using this analogy, the full system, containing multiple molecules loosely 
bound to each other near the wall P(7i, 72), may be thought of as a molecular cluster. While 
the black hole molecule may itself be made of atoms (i.e. single centered black holes), the 
nature of these constituents is irrelevant for the problem at hand. 

Our aim is to provide a wall-crossing formula which expresses Q~ in terms of Q + . In 
supergravity the difference Q~ — Q + is accounted for by the bound states of multiple black 
hole molecules carrying charges of the form Mtfi + N^, since they contribute to Q~ but not 
to Q + . Our goal in the rest of this section and the next section will be to use this viewpoint 
to derive the wall-crossing formula. 



r: {M 7 i + iV72, M,N>0, (M, N) ^ (0, 0)} , 



(2.5) 



(M71 + iV 72 , M'71 + N' l2 ) Im (Z Mjl+Nj2 Z M , J1+N , 72 ) < . 



(2.6) 
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2.3 Bose/Fermi statistics to Maxwell-Boltzmann statistics 

Let 70 be a primitive point on the charge lattice and let d(s7o) be the number of bosonic 
states minus the number of fermionic states of a black hole molecule 3 carrying charge S70 
moving in some appropriate background. We shall consider a system carrying total charge 
tc-fo consisting of a gas of mutually non-interacting black hole molecules carrying charges 570 
for different integers s. A typical configuration will contain m s black hole molecules with 
charge 570, subject to the constraint 

^sm s = fc. (2.7) 

s 

The net contribution to the index from such configurations is given by 

1 (d(s7o) +m s - 1)!' 



e n 



{ m s} 
2s sm s =fc 



m s \ (d(s l0 ) - 1))! 



(2.8) 



For bosons d(s7o) > 0, and the above result follows from the fact that m identical bosons 
occupying d states produce a degeneracy of ds = d(d + 1) ■ ■ ■ (d + m — 1) /ml. For fermions 
d(sjo) < 0, and the result follows from the fact that m fermions occupying \d\ states have 
total degeneracy dp = (\d\)(\d\ — 1) • • • —m+l)/m\ and index (— l) m dF = d(d+l) ■ ■ ■ (d + 
m — l)/m\. It will be convenient to encode (2.8) into a generating function 

1 (d(s l0 ) + m s - 1)\ 



m s \ (d(s l0 )-l))\ 



_ x sm s 



= JJ(1 - x s )- d ^ . (2.9) 

s 

We shall now prove that exactly the same contribution to the index is obtained if for 
each 7 we replace ^(7) by 

J( 7 ) = ^m^dij/m) (2.10) 

m|7 

and treat the particles as obeying Maxwell-Boltzmann statistics rather than Bose or Fermi 
statistics. For this we calculate the generating function of the index of a gas of non-interacting 
Boltzmann black hole molecules carrying charges S70 for different integers s. This is given 
by 



9{x) = EII^ (d(s l0 )) ms x sm ° = exp HT d(s l0 )i 

{m s } s s ' \ s 

= exp J <i(s7o/m)m~ 1 x s I = exp I d(kj )m~ 1 : 

\ s m\s I \ k m 

= exp ( - d ( k ^o) hi(l - x k ) J = JJ(1 - x k )- d{ ^ 0) . (2.11) 

V k J k 

3 We exclude from this counting the fermionic zero modes associated with broken supcrsymmctry. A 
formal way of doing this is to use helicity supertraces [44, 45] instead of the Witten index. 
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The m s \ in the denominator in the first line is the symmetry factor required for resolving 
the Gibbs paradox. Comparing (2.11) with (2.9) we see that the generating functions f(x) 
and g(x) are the same. Thus we are allowed to treat the black hole molecules as Boltzmann 
particles as long as we use the effective index ^(7). 

In general d(s'yo) receives contribution from the intrinsic index f2 + (s7o) of the black hole 
molecules and from the orbital degeneracy describing its motion in the background of other 
charges. The contribution to the orbital part of a black hole molecule of charge S70 is ex- 
pected not to be affected by the presence of the other black hole molecules carrying charges 
Z70 for any integer / since the symplectic product (kj , Ijo) vanishes and as a consequence 
the particles are mutually noninteracting. In particular the repulsive electromagnetic inter- 
actions cancel against the attractive gravitational and scalar exchange interactions. 4 On the 
other hand the orbital degeneracy is expected to depend on the background produced by 
other black hole molecules carrying charge not proportional to 70. We shall not need the 
explicit form of this contribution which is hard to compute in general when there are multiple 
other black hole molecules distributed in space, but use the fact that for a charged particle in 
a magnetic field the number of states per unit area in the lowest Landau level is proportional 
to the charge carried by the particle. To be more specific, we compare a configuration of 
a molecule of charge S70 moving in the background of other molecules carrying arbitrary 
charges with a configuration of s closeby molecules each of charge 70 moving in the same 
background. In this case the phase space volume element for the molecule of charge 570 will 
be s times the phase space volume element for each molecule of charge 70. Thus for a fixed 
background the orbital density of states for a black hole molecule carrying charge 570, being 
proportional to the phase space volume element, will be s times the orbital density of states 
of a molecule of charge 70. Thus we have g?(s7o) / '^(70) = s ^ + (s7o)/^ + (7o), an d hence 5 

d(s l0 )/d( l0 ) = ^m- 1 rf(s 7 o/m)/ C /(7o) = s ^ m~ 2 n + (s lo /m)/n + ( lo ) . (2.12) 

m\s m\s 

Comparing this with 

d( sl0 )/d( l0 ) = s n + (s l0 )/n + ( l0 ) (2.13) 

we see that replacing d(j) by ^(7) is equivalent to replacing fi + (7) by 

= m- 2 VL + (ci/m) . (2.14) 

We shall see in §4 that the fractional DT invariants 0(7) arise naturally in the KS wall- 
crossing formula. 

We end this section with a word of caution. For a generic interacting system of bosons 
or fermions the effect of statistics and interaction cannot always be disentangled to map the 

4 At short distance they may interact via the exchange of massive string modes and also via dipole and 
higher order interactions due to massless particle exchange, but we do not expect these interactions to affect 
the analysis of supersymmetric index. 

5 Note that in this argument the sign of d, and hence the statistics of the particle, is determined by the 
sign of il + . Although the orbital angular momentum contributes to the index of the final configuration, they 
arise from the electromagnetic field, and hence do not affect the statistics of the individual particles in the 
halo. 
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problem to an equivalent problem with Boltzmann particles. Consider for example the case 
of an interacting system of two identical particles for which at a certain energy eigenvalue E 
we have n s symmetric and n a anti-symmetric wave-functions. In this case we shall not get 
any simple map from the bosonic or fermionic system to a system of Boltzmann particles. 
Only if the identical particles are non-interacting so that multi-particle wave-functions can 
be constructed from (anti-)symmetric products of single particle wave functions, we have a 
simple map from a Bose/Fermi gas to a Boltzmann gas. 

2.4 General wall-crossing formula and charge conservation 

The analysis of §2.3 leads to the following prescription for computing wall-crossing from 
supergravity black holes. Suppose in the chamber c_ we have a configuration of multi- 
centered black hole molecules, consisting of m r s centers of charge (rji + S72) for different 
values of (r,s). These molecules interact via long range electromagnetic, gravitational and 
other massless particle exchange interactions. We now consider a quantum mechanics of this 
system regarding the different centers as distinguishable particles, each with unit degeneracy, 
and denote by §({m r]S }) trace of (— l) 2Js in this quantum mechanics. Then the wall-crossing 
formula will be given by 

Afi(M7i + N^ 2 ) = 1T(M7i + N>y 2 ) - ^(M^ + N^ 2 ) 
Yl 9({m r ,s})Yl 



{>™r,s} r,s 



(2.15) 



For gcd(M, N) = 1 we have ^(M^ + N^) = ^(M^+N^)- Eq. (2.15) then follows from 
the fact that the left hand side represents the change in the index and the right hand side 
represents the total contribution from the bound states of black hole molecules which exist in 
the chamber c_ but not in the chamber c+. For gcd(M, N) > 1 the indices (M71 + N~/ 2 ) 
are somewhat ill defined since the total index in the sector of charge (M71 + A^) can 
receive contribution also from unbound multiparticle states carrying charges (M71 +N r y 2 )/ s 
for integers s \ (M,N). Thus the unambiguous quantity is the effective index which enters 
the formula for the index of a bigger system of which the system with charge (M71 + A72) 
may form a subsystem [22]. This is what we call f2 ± and this is the quantity whose jump is 
computed by the right hand side of (2.15). 

A slightly different way of expressing (2.15) is: 

«-w-^w = E E.pd&IC^- (2 ' 16) 

n ^ 2 {ai,...,a„ef } 
7=«iH hOn 

where the sum runs over all possible unordered decompositions of a into vectors a±, . . . , a n G 
T. The integer |Aut({aj})| is defined as follows. If all the a^'s are distinct then |Aut({«j})| = 
1. If on the other hand the set (ai, . . . , a n ) consists of mi copies of a vector m 2 copies 
of a vector f3 2 etc. then |Aut({aj})| = n a m a'- Finally g(cv 1 , . . . ,a n ) represents the index 
associated with an n-centered black hole configuration in supergravity, with the individual 
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centers carrying charges cti,a2, • • • , oi n and treated as distinguishable particles carrying no 
internal degeneracy. Note that by an abuse of notation we have used for the argument of 
g two different representations of the index of multi-centered black holes - one where the 
arguments are charges carried by individual centers and the other where the arguments are 
integers specifying how many of the centers carry a given charge vector rji + 572- 

An immediate consequence of (2.16) is 'charge conservation' - the sum of the charges 
appearing in the arguments of Q + on the right hand side of the equation is equal to the 
argument of AQ on the left hand side of this equation. In contrast if we had written the 
wall-crossing formula using the indices Q + on the right hand side then there is no manifest 
charge conservation. This is a consequence of the fact that the use of Q allows us to use 
Maxwell-Boltzmann statistics for computing the contribution to the index due to multiple 
black hole molecules. In contrast if we had used Bose or Fermi statistics then manifest 
charge conservation is spoiled by the symmetrization effect since the degeneracy of k identical 
particles carrying index fl + not only contains a term proportional to (Q + ) k but also other 
terms containing lower powers of Q + . 

2.5 Semi-primitive wall-crossing from Boltzmann gas of black hole molecules 

In this section we shall derive the semi-primitive wall-crossing formula by counting the index 
of a gas of black hole molecules carrying charges S72 for different integers s, forming a halo 
around another black hole molecule of primitive charge 71. We denote by /C72 the total 
charge carried by the black hole gas. As noted in section §2.3, for this calculation we can 
regard the gas as one obeying Maxwell-Boltzmann statistics as long as we replace the index 
f2 + (s7 2 ) of a black hole molecule carrying charge 572 by f2 + (s7 2 ). The orbital motion of a 
black hole molecule of charge 572 around a molecule of charge 71 produces states carrying 
angular momentum (| (71,572)! — l)/2, and gives a contribution (— l)^ 7l ' S72 ^ +1 | (71, 572)! to 
the index [9]. Taking into account this additional factor we get the total contribution to the 
index from a single black hole molecule of charge 372 to be 

( _ 1) <7 1 -7 2 >+i| (7l)S72 )|0+( S72 ). (2.17) 

Since we have already chosen (71, 72) to be negative we can drop the absolute value sign and 
the +1 from the exponent. Thus if the halo consists of m s black hole molecules of charge 
S72 then the net contribution to the index is 



(2.18) 



Summing over all possible values of m s subject to the condition Y2 S sm s = ^5 an d mul- 
tiplying this by the index Q + ( , ji) of the black hole molecule of charge 71 we get a total 
contribution fl + ( r y 1 ) f2haio(7i> N), where 



fihato(7l,#)= E 



_ {"is} S 
Y1 B sm s =JV 



(2.19) 



This is the jump in the index due to a bound state of a black hole molecule of charge 71 
and a halo of black hole molecules carrying charges 572 for different integers s. In order to 
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calculate the total change in the index in the sector of charge 71 + across the wall of 
marginal stability "P( 7 i, 72), we need to sum over all possible bound states containing a core 
of charge 71 + Z72 and a halo of total charge (N — l)j2- Thus we have to sum over several 
terms of the form (2.19) with 7l replaced by 71 + Z72 and N replaced by (N — I) for different 
integers I. This gives 

N-l 

Afi( 7l + N l2 ) = ^ + (7i + h2) fihaio(7i, N-l) (2.20) 
1=0 

where we have used the primitivity of 71 + £72 to replace fi ± (7i + ^72) by fi ± (7i + £72) for 
< £ < N. This can be formalized in terms of the partition function 

00 

Z ± (l,g) = ^n ± ( 7l + iV72)g 7V . (2.21) 

N=0 

Then (2.20) can be written as 

Z-(l,q) = Z + (l,q)Z ba}o ( 7l ,q), (2.22) 

where 

Zhalo(7l, q) = 7^ £4alo( 7 l> N) q 
N>0 

We shall see later that this agrees with the KS and JS wall-crossing formula restricted to 
the semi-primitive case (eqs.(4.29), (4.30) and (5.26)). For N = 1 using 712 < 0, we recover 
the primitive wall-crossing formula (1.1). 

To recover the semi-primitive wall-crossing formula of [9] from (2.22), (2.23), we proceed 
as follows. First of all we note that the relation (2.14) can be inverted as 

fi W = E J2^ d ) > (2-24) 

where fj,(d) is the Mobius function (i.e. 1 if d is a product of an even number of distinct 
primes, —1 if d is a product of an odd number of primes, or otherwise). Using the identity 
[ldLi(l ~~ q d )^ d ^ d = e~ q we can now express (2.23) as 

^o(7i,9) = II i 1 - (-l) fc7 V) fc ' 712 ' n+(kl2) ■ (2.25) 

fe>0 

Eqs.(2.22), (2.25) give precisely the semi-primitive wall-crossing formula of [9]. 
2.6 Refined 'index' in supergravity 

Kontsevich and Soibelman also analyzed the wall crossing formula for the motivic DT invari- 
ants, which are quantized versions of the numerical DT invariants fi( 7 ; y). They enumerate 
the Betti numbers of the moduli space of BPS-states for given charge 7 in the weak string 
coupling regime, whereas f2( 7 ) equals the Euler characteristics of this moduli space, up to a 



00 



exp 2^g s (-i) 



(7l> s 72> 



(71, S7 2 )fi + (s7 2 



(2.23) 



K s=l 
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sign. Physically, the motivic DT invariants keep track of the angular momentum quantum 
numbers carried by the black hole at weak string coupling where the system may be repre- 
sented as a collection of D6-D4-D2-D0 branes wrapped on a Calabi-Yau 3-fold [33, 34]. A 
simple way to do this is to introduce an extra factor of y 2j3 inside the trace in (2.4) besides 
the (— 1) 2J3 factor that is already present in this definition of the index. Thus at a given 
point in moduli space, the refined 'index' 6 fi re f(7, y) is defined by [10, 33, 34] 

Orf(7,2/) = Tr^ {7) (-y) 2j3 = J^i-yT^M , (2-26) 

where Tr' denotes the trace over BPS states computed after removing the contribution 
from the fermion zero modes associated with the supersymmetries broken by the black hole. 
Alternatively we could compute the ordinary trace over all the BPS states and then divide 
the result by (2 — y — y^ 1 ) which represents the contribution from the fermion zero modes. 
The usual generalized DT invariants are obtained by setting y — 1, 

fi( 7 ) = ^(-l)"fi ref , n ( T ) = n iei ( 7) y = l) . (2.27) 

neZ 

In string theory (2.26) is not an index since it is not protected against quantum correc- 
tions as we switch on the string coupling. In supersymmetric gauge theories one can define 
an alternative version of this index as [21], 

OUT,!/) = Tr //(7) (2J 3 )(-l) 2J3 (-2/) 2/3+2J3 = Y,(-y) nn Un(l) (2-28) 

n&L 

where I3 is the third component of the SU (2)# symmetry. This is protected against quantum 
corrections. We shall however proceed with the definition (2.26) since our main interest is 
in string theory. Even though there is no general argument that protects this 'index' from 
changing as we vary the string coupling, and hence the DT invariants and the black hole 
degeneracies may not be equal, we may nevertheless expect that the structure of the wall- 
crossing formula at fixed coupling will remain the same. Thus we can regard the motivic KS 
formula as giving the change in f2 re f (7, y) across a wall of marginal stability at fixed value 
of the string coupling. With this in mind we shall analyze the jump in the motivic 'index' 
in supergravity and then compare this with the KS formula. Our supergravity analysis, 
compared to that in §2.3, will be somewhat heuristic; however the final result of this analysis 
will turn out to be consistent with the KS motivic wall-crossing formula. 

We begin by introducing some notations. We have already denoted by fi re f(7, y) the 
refined 'index' computed by introducing a weight factor of y 2j3 into the trace in (2.4). More 
generally we shall use the subscript re f to denote various quantities in which the trace over 
different angular momentum states has been performed with a weight factor of y 2J:i . With 
this the analog of (2.10) takes the form 

d{y,Js)= Yl m- 1 d(7/m,J 3 /m). (2.29) 

m\f,2J 3 

6 We shall use the words motivic and refined interchangeably. 
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A word of caution is warranted here. Since the full system is rotationally invariant, the 
states of this system can be characterized by their angular momentum. However when 
we examine the motion of one subsystem in the background of the other, the background 
generically breaks rotational invariance and hence states can no longer be classified by their 
angular momentum unless the background is generated by a point source (or a spherically 
symmetric source). We shall nevertheless proceed as if each subsystem consisting of a set of 
identical particles moved in the background produced by a point source so that an assignment 
of angular momentum quantum numbers to such individual subsystems were possible. Based 
on this assumption we shall arrive at an expression for the motivic index of the whole system 
in terms of the index carried by the individual molecules. This procedure can be justified a 
posteriori by the fact that it allows for a physical understanding of the motivic generalization 
of the KS wall-crossing formula. 

After multiplying (2.29) by y 2Ji and summing over J3 we get 

d«f(7,y)=X) E rn- 1 d( 1 /m,J 3 /m)y 2J ^^T, m ~ ld ^/ m ^ J ^y 2mJ ' 

J3 m\-y,2J3 m|7 Jg 

= J]m-V ref (7/m,y ro ). (2.30) 

m|7 

Our next task is to find the generalization of (2.12). Let us denote by c/ or b(7, ^3) the degen- 
eracy due to orbital motion of a black hole molecule of charge 7 in some fixed background. 
Again we pretend that the background is spherically symmetric so that it makes sense to 
assign definite angular momentum quantum numbers to the orbital states of individual sub- 
systems. Then we have 

drcf (7, V) = 0£f (7, 2/)4rb;rcf (7, 2/) , (2-31) 

where 

<Wf (7, </) = E d - b ^' Js ) y 2Jz ■ ( 2 - 32 ) 

Eq. (2.30) and (2.31) now give 

<Lf(7, y) = E m ~ l ^-f (7/m, y m ) d orh , tc{ ( 7 /m, y m ) . (2.33) 

We shall now try to express d ov ^ rc {{pf /m, y m ) in terms of G? or b ;re f(7, v)- F° r this (still pre- 
tending that we have a rotationally invariant subsystem) we shall decompose the orbital 
spectrum into SU (2) representations and denote by 6(7, J) the coefficient of the character 
of the representation of spin J . Then we have 

7/ 2J+l _ ,,-2.7-1 

<We f (7, v) = E J ) (y 2 J + y 2J ~ 2 + ■■■+ ^ 2J ) = E J ) — — • ( 2 - 34 ) 

j j y-y 

We now use the fact that for any m G Z + , and for spherically symmetric background, we 
have 

6(7, J) = 6(7717, J'), 2J' + 1 = m(2J+l). (2.35) 
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Effectively (2.35) follows from the fact that increasing the charge of the molecule by a factor 
of m changes the angular momentum carried by the lowest Landau level such that the 
degeneracy of the Landau level gets scaled by a factor of m. Using this we get 

; y(2J' + l) _ y~(2J' + l) y _ y-l 

dorb;r e f(l/m, y™) = ^ 6(7, J') ^ = 7^,. 17^ rf orb;ref (7, V) ■ (2-36) 



y'm y~ m 



y'm y~ m 



Substituting this into (2.33) we arrive at 

where the "rational motivic invariants" f2 re f are defined by 



(2.37) 



^rcf(7,y) = 



y - y 



m|7 



m (y m — y m ) 



Qref(l/m,y r 



(2.38) 



This shows that in computing the refined index of the full system we can treat the particles 
as obeying Maxwell-Boltzmann statistics provided we replace n re f by Q re f. As in the case of 
the classical DT invariants, the use of these invariants ensures that only charge preserving 
terms appear in any wall-crossing formula. 

The rational motivic invariants have appeared earlier in other contexts, for example in 
the construction of modular invariant partition functions in [29] . We would also like to point 
out that their structure is very similar to the free energy which arises in the computations 
by Gopakumar and Vafa [27]. The only difference is that in the Euclidean setup of [27], the 
factor (y d — y~ d ) on the right-hand side of (2.38) is replaced by (y d — y~ d ) 2 ■ As in this case, 
the generating function of the rational invariants fl re { leads to a product formula 



E 



y-y 



log 



n (i 



l+n+2£ -r<t>\ (- 1 )™ Q rof,n(7) 

y ) 



ngZ,«>0 
7 



(2.39) 



where is a vector of chemical potentials conjugate to the charge vector 7. Note that this 
product structure is lost in the limit y — > 1. 
The analog of (2.15) now takes the form 7 



A0 rcf (M 7l + N l2 , y) = ft- f (M71 + A^ 7 2, y) - ^tA M li + ^72, y) 

1 



E 

{"»r,s} 



0ref({nV,«},y) J 



(^(^71 + 572,2/))^ 



(2.40) 



where g re {({m r)S } , y) measures Tr (— y) 2J:i from orbital motion of a set of distinguishable 
particles, containing m r>s number of particles carrying charges rji + 572- Similarly the 



7 Since a single centered BPS black hole is expected to carry zero angular momcntum[46, 47, 48], one 
might naively expect f2^ f (7, y) to be independent of y. However as discussed in §2.2, we allow for centers 
which consist of multi-centered black holes whose relative separation remains finite as we approach the wall 
of marginal stability. As a result f2 ro f (7, y) can be a non-trivial function of y. 
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analog of (2.16) is 



«>2 {«!,.. .,a„er} 

7=aiH ha,i 

consistently with the charge conservation property of the motivic wall crossing formula when 
expressed in terms of the combinations O ref . Note that even though our derivation of (2.41) 
has been marred by unreasonable assumption of spherical symmetry in the dynamics of 
various subsystems, each term in (2.41) is defined unambiguously so that it can be put to 
test against known results. 

We can also easily derive the semi-primitive version of the motivic wall-crossing formula 
by following the logic of §2.5. It takes the form 8 

Afi ref ( Tl + iV7 2 ,2/) = 

^Z} _ _,,_-.[• 1 / (71, S72> _ (_ v )-(ni,sn)\ _ \ m H 

E^Wft + «•».»> E II ^ (( ' V) jK^y)) 

1=0 {m a } s 1 s ' v v y y / / j 

E s sm s =JV — I 

(2.42) 

We shall see later that this is in perfect agreement with the prediction of KS motivic wall- 
crossing formula (4.72). 



3. Multi-black hole bound states and quiver quantum mechanics 

In order to have a complete wall-crossing formula we need to find explicit expressions for the 
functions #({«;}), g TC f ({««}, y) appearing in eqs.(2.16) and (2.41) respectively. This requires 
finding the spectrum of supersymmetric bound states of multi-black hole configurations in 
supergravity. As argued by Denef[6], the supersymmetric quantum mechanics of multi- 
centered BPS configurations can be viewed as the "Coulomb phase" of a quiver matrix 
model, valid at strong coupling. At weak coupling, the wave function is instead supported 
on the Higgs branch. It should be possible to compute the function g({ai}), g xe { ({ctj}, y) 
from either description. In this section we shall describe both these approaches. As we shall 
see, the description on the Higgs branch, described in §3.1, allows us to solve the problem 
completely. On the other hand the description on the Coulomb branch, described in §3.2, also 
gives a complete algorithm for finding g({ai}), but it is more difficult to solve it explicitly. 
Furthermore with some guesswork we can also arrive at a specific proposal for g ve { ({«i}, y) 
from the analysis on the Coulomb branch. 

Without any loss of generality we can arrange the a^'s so that 

ccjj = (cui, ctj) > for i < j . (3.1) 

8 Notc that in the semi-primitive case our heuristic derivation becomes rigorous since identical particles 
carrying charge S72 move in the spherically symmetric background produced by the charge 71 + Z72 . Since 
[22] argues that general wall-crossing formula can be derived from the semi-primitive formula, we can use 
this to justify our general claim (2.41). 
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Now if we represent a vector Mji + in f 1 by the point (M, N) in the Cartesian co- 
ordinate system, then in this representation a pair of vectors (a, (5) will follow clockwise 
(anti-clockwise) order if (a, /3) is positive (negative). The condition (3.1) then implies that 
the vectors a±, ■ ■ ■ a n are arranged in a clockwise fashion. Throughout this section we shall 
work with this particular ordering of the a^s. 

3.1 Higgs branch analysis 

As has been argued by Denef [7] , the bound state spectrum of multi-centered black holes can 
also be computed using quiver quantum mechanics. For computing g rci ({c^}, y) we need to 
study the bound states of n distinguishable particles. In this case the quiver takes a simple 
form with n-nodes each carrying a U(l) factor, and a^- arrows from the node i to the node 
j for i < j. In particular since the arrows always go from lower to higher node, there are no 
oriented loops. 

Now for quivers without oriented loops, Reineke [39] has computed the generating func- 
tion of the Betti numbers of the space of semi-stable solutions to the D-flatness conditions. 
Physically they determine the number of supersymmetric quantum states carrying given an- 
gular momentum J3. The formula takes a simple form when all nodes carry U(l) factors 
and we shall state the formula for this special case. According to this formula Tr ((— y) 2j3 ), 
which can be identified with the function g re f({aj}, y), is given by 

g rc ^i,...,<y n ,y) = (-y)- L (y 2 -l) 1 - n (_i)-y E^E^^^'-f ? (3 . 2) 

partitions 

where the sum runs over all ordered partitions of (ai + - ■ - + a n ) into s vectors 0( a > (1 < a < s, 
1 < s < n) such that 

1. £a0 (o) =«! + •" + °n- 

2. /?( a ) = Yli m i a i wr th m!f^ = or 1 for each a,i. 

3. arg \ Y^ b a=1 Z^a)^ > arg (Z ai+ ... +an ) for all b between 1 and s — 1 in the chamber c~. 

Using the fact that in the chamber c~ arg(Z a ) > arg(Z (3 ) implies (a, (3) > 0, we can express 
condition 3. as 

{a \ a i + ■■■ + a nj > v b with 1<6<s-1. (3.3) 
In (3.2) L is a constant given by 

L = (l-rO + X>i, (3-4) 

i<j 

in such a way that (3.2) is invariant under y — > y~ l . Physically L represents the maximum 
2J 3 eigenvalue that the system can carry. 



- 19 - 



We shall now illustrate how this formula works by computing g re f ({««}, y) for n = 2,3 
and 4. First consider the n = 2 case with a± 2 > 0. In this case we have two possible ordered 
partitions satisfying the conditions stated above: 

{a 1 + a 2 }, {011,012}. (3.5) 

The first term contributes y 2ai2 and the second term contributes —1 to the sum. In total, 

g Icl K a 2 , y) = (-y) 1 ^ (y 2 - l)' 1 (y 2a - - 1) = V — - y ~™ . (3.6) 

y-y 

This correctly produces the spectrum of two centered bound states. Taking the y — > 1 limit 
gives 

g(a u a 2 ) = (-l) Ql2+1 « 12 . (3.7) 

Let us now turn to the n = 3 case. We take ai, a 2 , a 3 to be three distinct elements 
(not necessarily primitive) of Y such that a\ 2 , aii 3 and a 23 are all positive. For definiteness 
we shall choose the oij's such that a>i 2 > oc 23 . In the convention described above (3.1), the 
charges listed below follow a clockwise order as we move from left to right: 

ai, (cti + a 2 , a 1 + a 3 ), oli + a 2 + a 3 , a 2 , a 2 + a 3 , a 3 . (3.8) 

The relative ordering of the vectors inside ( ) is not determined unambiguously but is unim- 
portant. The condition (3.3) is now easy to implement: for every b, Yl a =i must be one 
of the vectors to the left of a\ + a 2 + a 3 in the list (3.8). In this case the possible ordered 
partitions of ol\ + a 2 + a 3 satisfying (3.3) are: 

{«i+a 2 +« 3 }, {ai, 0:2+0:3}, {®i+0! 2 ,a 3 }, {«i+a 3 ,a 2 }, {«i, o 2 ,o 3 }, {ai,a 3 ,a 2 }. 

(3.9) 

This gives, after a cancelation between the second and the last contribution, 

g Iei (ai, a 2 , a 3 , y) =(-i)°>»+°™+°™ ( y - 

I '„,C«12+ai3+Q!23 „,ai2— "23— "13 „,ai3+«23 — "12 I ,.— Q!12— «13— "23 )^ 

\y -y -y +y ) (3.10) 

= (-l) ai2+Ql3+Q23 sinh(z/a 12 ) sinh(i/(ai 3 + a 23 )) , 

sinn v 

where v = Iny. Taking the y — > 1 limit we get 

g(a u a 2 , a 3 ) = (-l)^+^+a 23 ^ ^ + ^ (3 n) 

If instead a-i 2 < a 23 , a similar reasoning leads to 

Steffi, 012, a 3 , y) = (_i)°i2+ai3+"23 . 1 s inh(^a 23 ) sinh(z/(ai 2 + ct 13 )) , 

sinn v 

^ (-1)^+^3+^23 ^ (ai2 + fti3) (3 12) 

Next we consider the case n = 4. We choose 4 vectors a^, a 2 , a 3 , a- 4 such that in the 
convention described below (3.1), different linear combinations of the a-j's follow the following 



- 20 - 



clockwise order as we move from left to right in the list: 



CKi , (ai + a 2 , Oil + Ol 3 , Oil + Ol 2 + Ol 3 ) , Ol 2 , («2 + "3 , "1 + "2 + « 4 ) , 

011 + oi 2 + a 3 + "4 , "35 «i + a 3 + a 4 , (ai + a 4 , a 2 + a 3 + a 4 , a 2 + a 4 , 03 + "4) , "4 , 

(3.13) 

where again, the order of terms between brackets is irrelevant. We list below the allowed 
partitions consistent with the three conditions described above and the corresponding con- 
tribution to the summand: 



{cm + oi 2 + a 3 + a 4 } 


^2(012+013+014+023+024^ 


-"34) 


{ai,a 2 + oi 3 + on} 


_y2(o 2 3+0 2 4^ 


-034) 


{a 2 , oil + a 3 + a 4 } 




-034) 


{«i + a 2 , oi 3 + a 4 } 




-034) 


{aii + a 3 , a 2 + a 4 } 


_y2(«i3+a23" 


-024) 


{«! + a 2 + 0:3, "4} 


_y2(«i2+ai3- 


-Q!23) 


{a 2 + a 3 , ai + a 4 } 


_y2(ai2+ai3+ai4- 


-023) 


{«! + a 2 + a 4 , a 3 } 


_y2(ai2+«i4+a 2 4- 


-O34) 


{ai,a 2 ,a 3 + a 4 } 


^2a 34 


{«1, «2 + «3, «4} 


^2d23 


{«i + « 3 , a 2 , a 4 } 


y2(ai3^ 


-«23) 


{«i,a 3 ,a 2 + a 4 } 


^2(a 2 3- 


-O24) 


{a 2 ,a 3 ,aii + a A } 


^2(ai2+ai3^ 


-Q14) 



oi 2 + a 4 , a 3 } 


^2(024+034) 


{a 2 , ai, a 3 + a 4 } 


^2(^12+034) 


{a 2 , «i + a 3 , a 4 } 


^2(012+013) 


{a 2 , «i + a*, a 3 } 


^2(012+014+034) 


{ai + a 2 , a 3 , a 4 } 


^2ai2 


{ai + a 2 , a 4 , a 3 } 


^2(ai2+«34) 


{a 2 + a 3 , ai, a 4 } 


y2(«i2+ai3+«23) 


{«i, a 2 , a 3 , a 4 } 


-1 


{«i, a 2 , a 4 , a 3 } 


_^2a34 


{ai, 03, a 2 , a 4 } 


^2Q23 


{a 2 , ai, a 3 , a 4 } 


^2qi2 


{a 2 , 03, ai, a 4 } 


_^2(Qi2+Ol3) 


{a 2 ,«i,a 4 ,a 3 } 


_y2(ai2+«34) 



Adding these terms and substituting into eq.(3.2) we arrive at 



g re {(ai,a 2 ,a 3 ,a4,y) 

_( -|\012+013+014+023+024+034+1 



(y-y^ 1 ) 3 

012+013+014—023—024 — 034 „,Ol2 +013 +014 +023 — 024~ 034 



x [y 

«i2+ai3+ai4— 023—024+034 



y 



+ y 



012—013+014—023—024+034 



v 
v 



-O12— Q13— O14 — Q23 — Q24— 034 



-O12+Q13— 014+Q23+024— O34 



^Ol2-Oi3+Oi4— 023+024+034 _|_ ^Oi2+Oi3+Oi4+Q!23+024+034 
_|_ ^—012+013— ai4 +023—024—034 _|_ ^—012—013 — 014+023+024—034 

_|_ ^—012—013—014—023+024+034 ^ — 012 — 013 — 014 +023 +024 +034^ 



\Ol2+Ol3+Ol4 +023 +024+034+l 



1 



sinh v 



sinh(z/«i 3 ) sinh(^(-a i2 + a 23 + a 24 )) sinh(z/(a 14 + a 34 )) 
+ sinh(i/«i 4 ) sinh(i/a 23 ) sinh(z/(-o+ 2 - a 13 + a 24 + a 34 )) 
+ sinh(^ai 2 ) sinh(^(a 13 + a 23 )) sinh(^(ai 4 + a 24 + a 34 )) 



(3.14) 
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Taking the y — > 1 limit we get 



g( y a 1 ,a 2 ,a 3 ,a 4 ) = (-l) 1+ ^<,^ x 

["l2 "l3 "24 + "l3 "l4 "24 + "l2 "23 "24 + "l4 "23 "24 

+"12 "23 "34 + "13 "23 "34 + "14 "23 "34 + "l3 "24 "34 ] • (3. 15) 

We have also derived the analog of (3.14) for n = 5 but we suppress the result for the sake of 
brevity. Similar results can be obtained for different choices of orderings, while cases where 
some of the final states coincide can be obtained by taking suitable limits. For example 
consider the case g re f ("1, "2, "2, y) with ayi > 0. This can be considered as a special case of 
(3.10) in the limit a 3 — > a 2 and gives 



We shall later verify that the various explicit results given in this section are in perfect 
agreement with both KS and JS wall-crossing formulas. 

Thus, eq. (3.2) provides a complete algorithm for computing the coefficient g re f ({«i}, y) 
in (2.41) for any number of ttj's. This result is based on the study of multi-centered bound 
states in supergravity, even though we had to rely on mathematical results about moduli 
spaces of quiver representations. The key point is that we only needed invariants of Abelian 
quivers to compute the dynamics of internal (or hair) degrees of freedom, while each center 
could still be regarded as a macroscopic solution (provided its charge is large enough). In the 
next subsection, we relinquish the gauge theoretical crutch provided by Reineke's formula 
and directly quantize the internal degrees of freedom of the bound state. 

3.2 Coulomb branch analysis 

In this subsection we shall try to reproduce the results of §3.1 by directly quantizing a con- 
figuration of multi-centered black holes. We begin by reviewing some relevant properties of 
these solutions. Consider a supersymmetric solution describing n black holes, with different 
centers carrying charges ai, . . . , a n £ F located at fi, . . . , r n . We shall define 



fiW("l,"2,"2) 



1 



sinh(i/a 12 ) sinh(2z/a 12 ) . 



(3.16) 



1 ij — I ' i 1 j I • 



(3.17) 



The equations governing the locations fi are [6] 




(3.18) 



where 



c 4 = 2Im(e-%) , 



(p = arg(Z ai+ ... +Q J . 



(3.19) 
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Here Z a denotes the central charge for the charge a, computed with the asymptotic val- 
ues of the moduli fields. The constants q depend on the moduli through Z a) and satisfy 
V c- = 9 

First we shall show that in the chamber c + none of these solutions exist. For this note 
that at the wall of marginal stability ^(71,72) the central charges of 71 and 72, and hence 
of all the vectors a«, become aligned. As a result, near this wall the real parts of e~^Z ai are 
all positive; we shall denote these by A4. On the other hand it follows from (3.1) that for 
2 = 1 the left hand side of (3.18) is positive and for i = n it is negative. Thus we must have 

e-**Z ai = At + iBx, e-*+Z an =A n -iB n , A 1 ,B 1 ,A n ,B n >0. (3.20) 

This gives a necessary condition for the multi-centered solution to exist, 

(ai, a n ) Im (Z ai Z an ) = (a u a n ) (A 1 B n + A n B x ) > . (3.21) 

On the other hand the chamber c + has been defined such that the right hand side of (3.21) 
is negative (see (2.6)). This shows that a multi-centered solution of the type analyzed above 
cannot exist in the chamber c+. Note that this also proves that scaling solutions [7, 49, 9], 
whose existence does not depend on the moduli, cannot exist for charges a\, . . . , a n G T. 10 

From now on we work in the chamber c_. For an n-centered configuration, the location 
of the centers is specified by 3n coordinates fi (1 < i < n). Removing the trivial center 
of mass degrees of freedom we are left with 3n — 3 coordinates. For multi-centered BPS 
solutions, the relative distances = |ff — fj\ must satisfy (3.18) for i = 1, . . . , n. These 
equations are linearly dependent, since the sum over % is trivially satisfied. This gives (n — 1) 
independent constraints. The moduli space of multi-centered solutions is then a (2n — 2) 
dimensional space A4 n (ci, . . . ,c n ). In the case of interest here, where all ccj lie in T, M. n is 
compact. 

As shown in [40], for fixed values of Cj, Ai n carries a symplectic form u given by the 
restriction of the two-form 11 

Kj J 

from ]R 3n_3 to the moduli space -M n (ci, . . . , c n ). This symplectic form is invariant under 
SU(2) rotations. The moment map associated to an infinitesimal rotation 5r = e* A r*is just 

9 In the original analysis of [6] each center was regarded as a single centered black hole. We shall use a 
slightly more general interpretation in which each center is allowed to be a black hole molecule - composite 
of two or more single centered black holes with charges of each center lying outside T. By working at a 
point close to the wall 7 , (7i,72) we can ensure that the distance between these molecules, denoted by r,-j in 
eq.(3.18), is much larger than the internal size of each molecule, and hence (3.18) is a valid description of 
the inter-molecular distance for this configuration. 

10 We can also directly see this as follows. For the scaling solutions the right hand side of (3.18) vanishes [9]. 
Now for i = 1 all the terms in the left hand side of (3.18) are manifestly positive due to the choice (3.1), 
making it impossible to satisfy this equation. 

n Our normalisation differs from a factor of two from the one used in [40]. This ensures that lo/(2tt) has 
integer periods. 
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e • J, where 



IE^ct ( 3 - 23 ) 



7 2 ^ " w |r 

is the angular momentum. Thus the spectrum of supersymmetric bound states can in prin- 
ciple be obtained from geometric quantization of this phase space and the information on 
angular momentum, required for computing g Te f, can be found by studying the J 3 eigenvalues 
of these bound states. 

We now review the results of [40] for the bound state spectrum of 3-centered black holes 
in the chamber c_. 12 In this case we have (aii, ay) Im (Z ai Z aj ) > and hence the clockwise 
ordering of the ajj's will imply clockwise ordering of the Z a /s. Furthermore we shall restrict 
the aij's to satisfy 

"12 > "23 (3.24) 

so that the clockwise ordering of the a^'s and their various linear combinations are given by 
eq.(3.8). Thus the same is true for the corresponding Z's. We now explicitly write out the 
equations (3.18) as follows: 

^ + ^=ci, -^ + ^ = c 2 , ^ + ^ = - C3 , c^2Im(e-^). (3.25) 
ri2 ri3 r 12 r 23 r 13 r 23 

Since e 1 ^ = Z ai+a2+a . 3 /\Z ai+a2+a3 \, we have c\ + c 2 + c 3 = 0. Furthermore since according to 
(3.8) Z ai and Z ai+a3 precedes Z ai+a3+as in the clockwise ordering while Z a2 and Z a3 follow 
it, we have c\ > 0, c\ + c 3 > 0, c 2 < 0, c 3 < 0. Thus we can parametrize the Cj's as 

ci = a, c 2 = — a + 6, c 3 = — 6, a, 6>0, b<a. (3.26) 

We can express the general solution to (3.25) as [40] 

«12 «23 «13 /q „x 

ri2 = t r, ^23 = t , r i3 = — — r r , 3.27 

A — b A — a a + b — A 

for some constant real parameter A. The range of A is restricted by the positivity of each 
Tij and also the triangle inequality satisfied by the r^-'s. The positivity of the r^'s together 
with (3.26) give b < a < A < a + b. To study the consequences of the triangle inequality we 
express them as 



«23 «13 





A 


— a a 4- 


■6- A 


«23 




«13 


CK12 


A — a 


a 


+ 6- A 


A-6 


«13 




CK12 


«23 



a + b — A A — 6 A — a 



>0, 
>0, 

>0. (3.28) 



12 Note that if the three charge vectors a^s do not lie in a plane, then we can compute the bound state 
degeneracy using attractor flow trees [6], e.g. we could first approach a wall where the system splits into 
a pair of molecules, one with charge ot\ + a?, and another with charge 013 and so the index will be given 
by (— l) Ql3+Q23+1 (o: 1 3 + a23)0(ai + a.?)Sl(a.z) . Then we can approach another wall where the system with 
charge a\ +0:2 breaks up into a pair of molecules of charges a\ and 0:2 with index (— 1) Q12+I ai2^(ai)f2(a2). 
But when the three o^'s are in the same plane spanned by 71 and 72, they all move away from each other 
at a similar rate when we approach the wall of marginal stability "P(7i, 72) and we need to solve the 3-body 
bound state problem explicitly. Similarly if we have n centers with their charges lying in a plane then we 
need to explicitly solve the n-body bound state problem. 
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We need to find solutions to these inequalities in the range a < A < a + b. We begin near 
A = a + e for some small e. At this point r23 is large and the last of eqs.(3.28) is violated. 
As we increase A, at some value the last equality is saturated when 

r 2 3 = n 3 + r 12 . (3.29) 

It is easy to see that at this point the other two inequalities hold and hence above this bound 
the allowed range of A begins. This continues till one of the other inequalities fail to hold. It 
is easy to see that second inequality continues to hold for A < a + b but the first inequality 
is violated beyond some value of A close to a + b when 

r\3 = n2 + r 2 s • (3.30) 

The allowed range of angular momentum is given by the classical angular momentum 
carried by the system in the two extremes. At (3.29) the points (r* 2 , ^1,^3) he along a line and 
hence the angular momentum is given by (ai 3 + a 23 — «i2)/2. On the other hand at (3.30) the 
points (r*i, r 2 , r 3 ) lie along a line and we have total angular momentum («i 3 + a 23 + ai2)/2. 
Thus we have 

J- = (a 13 + a 23 - a 12 )/2, J + = (a 13 + a 23 + a l2 )/2 . (3.31) 

As was shown in [40] in quantum theory the upper limit J + is shifted to J + — 1 and states 
of all angular momentum between J_ and J + — 1 occur exactly once. This gives [40] 

J +Z} ? ,2J+1 _ ,,-2.7-1 
ff»f(«l,«*«s;!/) = (-l^+C+ow — x 

ti. y - y 

= (-l)«i3+« 2 3+c* 12 l -^- 2 {y 2J+ - y 2J - - y~ 2J - + y~ 2J+ ) 

(y-y ) 

= (-l) a «+a 23 +« 12 _J_ sinh ( u ( ai3 + « 23 )) S inh(i/«i 2 ) , (3.32) 
smh v 

in agreement with (3.10). 

We shall now generalize Eqs. (3.32) to an arbitrary number n of centers. For this we 
shall first simplify (3.18). Since we are interested in the situation where the a^'s lie in a two 
dimensional plane we have 

2 Im (e-'^Z ai ) = A (a h ^ aj) = A ^ a„ , (3.33) 

for some positive constant A. This constant can be removed by a rescaling of the variables 
r*j, but we shall keep it in our subsequent equations. This allows us to express (3.18) as: 



n 



1.1 



3=1 v j& 



A^aij. (3.34) 



Our strategy will be to relate Tr (— l) 2Js |/ 2j3 to an integral over the classical phase space 
M. n of solutions to eq. (3.34). Now, (— 1) 2J3 is a rapidly varying function on Ai n and does 
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not have a smooth classical limit. Our experience with the quantum theory for n = 2 and 
n = 3 nevertheless suggests that it takes the same value over all the quantum states and is 
given by (— \Y Jmax ~ n+1 , where J max — Y2i<j a ij/2 is the maximum allowed classical angular 
momentum. On the other hand, for y close to 1, y 2Js is a slowly varying function over the 
classical phase space and one expects that for large |ay|, its quantum expectation value is 
well approximated by integrating y 2Js over the classical phase space. Thus, we are led to 
consider 

£cia Ssi cai(W,?/)= \ ' 1( rrr / e 2v J * a/ 1 " 1 , v = lny . (3.35) 

[ATT) \n — ly. j Mn 

This formula should well approximate the refined index g Te { at large | a^- 1 and y close to 1 , but 
could in principle be corrected in the full quantum theory. Our experience with the two and 
three centered cases, as well as an explicit evaluation of the 4-centered case presented below, 
suggests that at y = 1 the classical phase space integral (3.35) in fact agrees with the exact 
quantum index g({cti}). In addition, the same integral (3.35), after a minor amendment to 
be described shortly, appears to agree with the exact refined index g re f({aj}) for all values 
of y. 

Now using the localization theorem of [50], we can express (3.35) as a sum over contri- 
butions from fixed points of the Hamiltonian vector field associated to the moment map J3, 
i.e. rotations along the z-axis. Fixed points are therefore multi-centered black hole configu- 
rations in which all centers are aligned along the z-axis, in an appropriate order consistent 
with (3.34). Furthermore, since all the relative distances between the centers are fixed by 
(3.34), the fixed points are isolated. Thus, fixed points of J3 are labelled by permutations a 
of 1, 2, . . . , n such that the centers are arranged in a given order along the z-axis, satisfying 
z c (i) < z a (j) if % < j. In this case the constraint (3.34) takes the form 

n 

a a{i)a{j) , (3.36) 

2=1 

which is equivalent to the extremization of the "superpotential" 

W = -- ^2 si g n (i ~ i)a*(i)a(j) lo S \ Z "V) ~ Mi) I _ A ^2 ai i Zi ■ ( 3 ' 37 ) 
At such a fixed point the third component of the classical angular momentum is given by 

i<j 

The localization formula of [50] now gives 

ga^doihy) = (-l) E ^ Q - +n " 1 (2z/) 1 ""^' s (a) <* a(i)aU) ? (3.39) 

IT 

where Yl' a denotes sum over only those permutations for which a solution to (3.36) exists, 
and s(a) is the sign of the Hessian of the matrix representing the action of 2v J3 on the 



E 

2=1 



a «r(Q«T(j) 
Z f{j) ~~ Mi) 



sign(j 
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tangent space of Ai n at the fixed point. To compute s{a) we make a convenient choice of 
coordinates on Ai n . Without any loss of generality we can choose r\ to be at the origin. At 
a fixed point of the action of J3, all the other points are then along the z-axis. We now note 
that to first order the relative distances between the centers remain unchanged if we displace 
each of the fi for i > 2 in the {x — y) plane. Thus these (2n — 2) coordinates provide us with 
a convenient parametrization of the moduli space of the solution near this fixed point. Let 
us denote them by (2 < i < n). The action of J 3 on these coordinates is simply an 

independent rotation in the (xj, yi) plane for each i. The Hessian of 2z/J 3 is given by (2z/)" _1 
up to a sign s{cr). To determine the sign we note that in the coordinate system {xi,yi} 
introduced above, J3 and u take the form: 

J 3 = ^^2 a <rd)°(j) ~ \ M ij( x i x i + ViVj) H > u = ^ M ijdxi A dyj H , (3.40) 

i<j 

where ■ ■ ■ denote higher order terms and 

Mu = a u -r— 13 + a ik - — - — - , Mij = -a ir — - — - for 1 ^ j, 1, j > 2 . 

\ z i\ k ^ k > 2 \ Zk Z{ \ \ z i Zi \ 

(3.41) 

It is worth noting that the matrix is also equal to the Hessian of the superpotential 
(3.37) with respect to the n — 1 variables Z2, ■ ■ ■ Z n , with z\ being set to zero. The sign s(er) 
of the Hessian associated with the fixed point is thus given by 

s(a) = signdetM. (3.42) 

Although the prescription (3.42) appears to treat Z\ on a different footing than the other 
Zi due to the gauge condition Z\ = 0, one could just as well have computed s(a) using a 
symmetric gauge condition ^=1 n Zi = ^- in deed, the same sign s(a) can be obtained as 
(the opposite of) the sign of the determinant of the Hessian of the superpotential 

A " 

W = W + - y^Zi (3.43) 
n ^— ' 

i=l 

with respect to all Zi, i — 1 . . . n and to the Lagrange multiplier A. To see this, note that the 
Hessian of W with respect to (A, Zi, Zi=2... n ) is given by 

/ 1/n A\ 
M — 1/n d 2 W/dzj W 1 (3.44) 
\A T W[ Mj 

where both A and W\ are {n—1) dimensional row matrices, with A = {1/n, ■ ■ ■ 1/n) and W\ = 
{d 2 W/dzidz2, ■ ■ ■ d 2 W / dz\dz n ). By adding the third to (n+l)'th rows to the second row and 
third to (n+l)-th columns to the second column and using the fact that Y17=i d 2 W/dzidzj = 
due to translation invariance, we can bring (3.44) to the form 



(3.45) 




-27- 



where denotes an (n — l)-dimensional row matrix with all entries 0. From this we see that 
det M = — det M. 

If there are more than one solution of (3.36) for a given permutation a then the right 
hand side of (3.42) should be replaced by a sum of sign det M over all solutions. Numerical 
evidence indicates however that there is at most one fixed point for a given permutation. 
Moreover, it suggests that s(a) can be expressed in terms of the permutation a through 

This is easily proven for the special critical points of W with sign(z CT (j) — z a ^) = sign(i — j) 
and sign(2 CT (j) — z a (j)) = sign(j — i). These represent solutions of (3.36) corresponding to 
permutations (12 . . . n) and (n . . . 21). Since these correspond to the global maximum and 
global minimum of J3, respectively, the matrix Mjj is positive definite or negative definite, 
respectively, leading to s(12...n) = 1, s(n . . . 21) = (— \) n ~ l . As we shall see in §3.3, the 
result (3.46) is required for consistency with the Higgs branch derivation presented in §3.1. 
This also suggests that if there are more than one fixed points for a given permutation their 
contributions should cancel pairwise leaving behind the contribution from or 1 fixed point. 

We do not expect the classical formula (3.35) and hence (3.39) to reproduce the full y 
dependent quantum answer for Tr (—y) 2J:i - after all the quantization of angular momentum 
is not visible classically. This is apparent from (3.39): while the terms inside the sum involve 
integral powers of y = e u and hence are compatible with charge quantization, the overall 
factor (2u) 1 ~ n does not respect charge quantization. Comparison with the exact results (3.2), 
(3.32) suggests a remedy 13 : replace the factor of {2v) 1 ~~ n by (2sinhz/) 1-n = (y — y -1 ) 1 "™. 
In the y — > 1 limit this will approach the classical result in accordance with the earlier 
observation that in this limit the classical and quantum results agree. Thus our proposal for 
the full quantum version of (3.39) is 

g rci ({a t },y) = (-1)^ ^ +n ~\y - y' 1 ) 1 '^ s{a) y^<> a ^ . (3.47) 

This reduces the problem of computing the function g Te t to identifying which of the permu- 
tations a are consistent with (3.36). This is a tedious but straightforward procedure. Below 
we give the results for n = 3 and n = 4 for the same order of various linear combinations of 
the ctiS as given in (3.8) and (3.13). 

For n = 3 the detailed analysis of the configuration space was carried out in [40] some 
relevant details of which were reviewed earlier in this section. Two of the four collinear 
configurations are given in (3.29) and two others are given by (3. 30). 14 This gives the 
following order of the centers along the z-axis and the value of s(a), 

{1, 2, 3; +}, {2, 1, 3; -}, {3, 1, 2; -}, {3, 2, 1; +} , (3.48) 

13 Note added in v2: in our subsequent work [58], we derive this multiplicative renormalisation from the 
Atiyah-Bott Lefschetz fixed point formula for the equivariant index of the Dirac operator on M n . 
14 For each of (3.29) and (3.30) we have two configurations related by z — > — z symmetry. 
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leading to 



X | ^"12+ai3+Q23 yQl3+"23-ai2 ^012-023-013 _|_ y-012-013-023 



J 



(3.49) 

in agreement with the result (3.32) following from exact quantization of the 3-centered 
system. 

Next we test (3.47) by working out the result for n — 4. There are 12 fixed points, whose 
orders and the associated s(a) are given by 

{1,2, 3, 4;+}, {3, 1,2,4;-}, {1,3, 4, 2;-}, {1,4,2,3;-}, {3, 1,4, 2;+}, {2, 3, 4,1;-}, 
{1,4,3, 2;+}, {2, 4, 1,3;-}, {3,2,4,1;+}, {2,4, 3,1;+}, {4,2, 1,3;+}, {4, 3, 2,1;-}. 

(3.50) 

Eq.(3.47) now gives 

n (n, n< n< n ni\ — ( 1 ^012+013+014+023+024+034 + ! („, „,- 1n \-3 

g Te i\otx, OLi, 0:3, a4,y) — ( — L) [y — y ) 

x ^012+013+014+023+024+034 ^012-013+014-023+024+034 ^Ol2+«13+ai4-023— «24+Q34 

_ yai2+ai3+ai4+a23— "24— O34 _|_ ^Oi2-Oi3+0! 1 4-a23-O24+O34 — Q12— Qil3— ai4+Q23+Qi24+a34 

_|_ ^Oi2+Oi3+Qi4 — 023— 024 — 034 y — "12+013— 014+023+024— 034 _|_ ^—012—013 — 014—023+024+034 

_|_ ^-012-013-014+023+024-034 _|_ ^-012+013 — 014+023—024—034 ^-012—013-014-023—024-034 

(3.51) 

in agreement with (3.14). 

3.3 Comparison of the results of Higgs branch and Coulomb branch analysis 

To compare the results of the Coulomb branch analysis described in §3.2 with the Higgs 
branch computation described in §3.1, note that the power y^<3 a "(i)"G) in (3.47) matches 
the power of y in (3.2), provided the ordered decomposition (3 is chosen as follows: break 
the sequence {a(i),i = 1, . . . , n} into increasing subsequences {cr(i a _i + 1), . . . , a(i a )}, = 
z'o < %\ < ■ • • < i s = n, where s is the number of such increasing subsequences. The vectors 
(3^ are then 

/3 (s+1 " a) = 2 ■ ( 3 - 52 ) 

i=i a -i+l 

For example for the permutation <r(123) = (132) there are two possible choices for the 
increasing subsequences: {{13}, {2}} and {{1}, {3}, {2}}, and the corresponding partitions 
are {02, ot\ + a 3 } and {02,03,01}. To see that the partition (3.52) generates the correct 
power of y associated with the permutation cr, we note that (3.52) implies that m s i +l ~ a = 
X)}=i„_i+i ^i,o-(j) * n (3-2). Leaving aside the factor of (y — y~ 1 )~ n+1 , the power of y in (3.2) 
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is now given by 



E^+ 2 E E E 



Ot<T(l)<T(k) 



E^+ 2 E 



a CT (0o-(fc) j 



(3.53) 



o-(i)<o-(fe) 



where in the second step we have used the fact that if k, I belong to different subsequences 
labeled by a and b then a > b will imply k > I. On the other hand if they belong to the 
same subsequence then, since the subsequence is increasing, the condition a{l) < a(k) will 
imply / < k. The right hand side of (3.53) gives precisely the power of y associated with the 
permutation a as given in (3.47). 

The sign associated with the partition (3.52) in (3.2) is (— l) s_1 . If the increasing 
subsequences are maximal, i.e. chosen so that it is not possible to build bigger increasing 
subsequences, then this sign is in agreement with the rule given in (3.46). For example, for 
the permutation cr(1234) = 3142 the maximal increasing subsequences are {{3}, {14}, {2}}. 
Thus it corresponds to the ordered decomposition {a 2 ,ai + ct 4 ,a 3 }, contributing with a 
positive sign. 

Now, there are typically many increasing subsequences associated with a given per- 
mutation, obtained by breaking up the maximal increasing subsequences into smaller in- 
creasing subsequences, and contributing with different signs. Thus, in order to determine 
if a given permutation contributes and with what sign, we have to combine the contri- 
butions from these different terms. In the previous example the other possible increasing 
subsequence is {{3}, {1}, {4}, {2}}, but the corresponding partition {a 2 , a 4 , a\, a-s} does not 
satisfy the condition (3.3) and hence does not contribute. In general however the situa- 
tion is more complicated. Consider for example the permutation <r(1234) = (3412). The 
corresponding increasing subsequences are {{34}, {12}}, {{3}, {4}, {12}}, {{34}, {1}, {2}}, 
and {{3}, {4}, {1}, {2}}, associated to the partitions {a.\ + a 2 , a 3 + a 4 }, {«i + a 2 , a 4 , a 3 }, 
{a 2 ,ai,a 3 + a 4 } and {a 2 , a±, a 4 , 03}, respectively. For the order given in (3.13) only the 
first and the third partitions are allowed by the rules (3.3). They contribute with opposite 
sign making the contribution vanish. This explains why the permutation (3412) is absent 
from the list (3.50). 

It is easy to convince oneself that all possible partitions of the vectors (a±, . . . ,a n ), 
whether or not they satisfy the condition (3.3), are in one to one correspondence with the 
set of all increasing subsequences of all the permutations of (12 . . . n) via the rule (3.52). So, 
the complete Higgs branch contribution can be generated by beginning with the maximal 
increasing subsequences associated with a given permutation and combining them with the 
contribution from other increasing subsequences associated with the same permutation. The 
following observations are now in order: 

1. If a given partition is not allowed by the rule (3.3) then all its subpartitions are also 
disallowed. For example in the example of the previous paragraph, once we know that 
{«i + a 2 ,a 4 ,Q!3} is not allowed, we can immediately conclude that {a 2 , a±, a 4 , a 3 } is 
also not allowed. 

2. If there are two or more maximal increasing subsequences of length two or more, 
then each of the maximal increasing subsequences can be independently broken up 
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into smaller increasing subsequences. The compatibility of a partition of a particular 
maximal increasing sequence with the condition (3.3) or not is independent on the 
partitioning of the other maximal increasing subsequences. For example for the par- 
tition {a i + a 2 ,a 3 + a 4 }, the compatibility of the splitting of a 3 + a 4 into {a 4 ,a 3 } 
can be determined independently of whether ol\ + a 2 is kept as a single element, or 
has been split into {ct2, ati}. For this reason we can associate, to each maximal in- 
creasing subsequence, a weight given by a sum of ±1 for each of the allowed splittings 
of that subsequence (+1 for splitting into odd number of subsequences, including the 
original maximal increasing subsequence, and —1 for splitting into even number of 
subsequences). The final weight is given by the product of the weights computed from 
each maximal increasing subsequence. 

In the example above the weight factor associated with a 3 + a 4 is 1 since it cannot be 
split, while the weight factor associated with a± + «2 vanishes since it allows a split 
{a 2 ,ai} with opposite sign. As a result the net weight is 1 x = 0. 

3. The problem of determining the contribution from a given permutation now reduces 
to computing the weight factor associated with each maximal increasing subsequence 
of that permutation. This can be done as follows. We begin with a particular maximal 
increasing subsequence and first consider all possible partitions of this subsequence into 
two smaller increasing subsequences. For this we need to simply insert a comma at one 
place that indicates how we divide the original subsequence. Not all such subsequences 
may generate partitions allowed by (3.3); let us assume that there are k possible places 
where we are allowed to insert the comma. This gives k terms, each with weight — 1. 
Now consider the possible partitions of the same maximal increasing subsequence into 
three increasing subsequences. It follows from the rule (3.3) that the allowed partitions 
are obtained by inserting a pair of commas into two of the same k possible positions. 
Thus there are kik — l)/2 possible terms, each with weight 1. This generalizes to 
partitioning into arbitrary number of increasing subsequences. Thus the net weight 
factor is 1 — A; + (2) — (3) + • ■ ■ = (1 — l) fc . This shows that the weight factor vanishes 
for k > 1, and is 1 for k = i.e. when it is not possible to subpartition an increasing 
subsequence satisfying (3.3). 

This leads to the following simple rule for deciding when a given permutation contributes 
and the sign of the contribution: A given permutation contributes if its maximal increasing 
subsequences generate a partition satisfying (3.3) via (3.52), and none of the other (non- 
maximal) increasing subsequences generate an allowed partition. The sign of the contribution 
is given by (3.46). 

We should of course keep in mind that the algorithm described above is not an indepen- 
dent result derived from the Coulomb branch, it is required to ensure that the Coulomb and 
the Higgs branch results agree. It will be interesting to find an independent derivation of 
this from the Coulomb branch analysis by directly examining the condition for the existence 
of solutions to (3.34). 

Before concluding this section we shall demonstrate how the algorithm given above can 
be used to give a derivation of the semi-primitive wall-crossing formula. We choose a n to be 
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71 and a\, ■ • ■ a n -\ to consist of m\ copies of 72, m 2 copies of 272 etc. in some fixed order 
(which can be decided by adding some arbitrary small vector to each of these charges which 
will be taken to zero at the end). Thus we have n — 1 + J2 S m s and J2i<j a ij = ~7i2 J2 S sm s- 
Since (71 + Ary 2 , li + Yls sm sl2) < for k < J2 S sm si ^ follows from the (3.3) that 71 must be 
part of the last partition. Consider now a permutation of 1, • • • n. In order that 71 is part of 
the last partition, the first maximal increasing subsequence in this permutation must contain 
the element n as its last element. This subsequence cannot be partitioned into smaller 
increasing subsequences since then ji will not be part of the last partition. Furthermore 
in order that the permutation gives a non-vanishing contribution the rest of the maximal 
increasing subsequences must each have length 1, since any maximal increasing subsequence 
of length 2 or more can be partitioned into smaller increasing subsequences without violating 
(3.3) and the result will vanish. This implies that the rest of the elements must be arranged 
in decreasing order in the permutation a. Thus the only freedom in choosing the permutation 
is in deciding which elements are part of the first increasing subsequence. Let there be k s 
copies of 572 in this set; this can be chosen in a total of Y[ s (™ s ) ways. The total number of 
partitions associated with this permutation is 1 + Y2s( m s ~ Thus the net contribution 
to g re{ is given by 



m ^ E * m * v - 1)- ^ n e r ) ( - i y 

Sj (3.54) 



\m s -k S y-2'y 1 2k s s 



n 



s k a =0 

(71,372} _ (_jA-<7i>s72>' 



in agreement with (2.42). 



4. Wall crossing from the Kontsevich-Soibelman formula 

The first solution to the problem of determining g({ai}) was given by Kontsevich and Soibel- 
man [10] and also independently by Joyce and Song [13, 14]. In this section we shall review 
the results of [10] and compare them with our results. Subsection 4.1 states the KS wall- 
crossing formula, Subsection 4.2 explains charge conservation and the following subsections 
apply the KS formula to determine the jump of ^(7) in various cases. In subsections 4.7 
and 4.8 we describe generalization of the KS formula to the motivic index. 

4.1 The KS formula 

To state the KS formula, we introduce the Lie algebra A spanned by abstract generators 
{e 7 ,7 G T}, satisfying the commutation rule 

[e 7 ,e y ] = K((7,7))e 7+7 /, (4.1) 

where we defined 

K (x) = (-l) x x . (4.2) 
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At a given point in moduli space labeled by the parameters {t a }, we introduce the operator 

r/ 7 (o = «p(n( 7 ;OE^) (4-3) 



d=l 



in the Lie group generated by A. The KS wall-crossing formula [10, 16] states that the 
product 

Km= n u ^ ( 4 - 4 ) 

7=^71+^72, 
M>0,N>0 

ordered so that as we move from the left to the right the corresponding Z 7 's are ordered 
clockwise, i.e. according to decreasing values of arg(Z 7 ), stays constant across the hyperplane 
of marginal stability V (71,72)- As t a crosses this locus, Q(7;t a ) jumps and the order of the 
factors is reversed, but the operator A 7l;72 stays constant. Thus, the KS formula may be 
stated as the equality 

n ^ i&7i+jv "72 = n u Mli+N12 , (4.5) 

M>0,N>0, M>0,N>0, 
M/Nl Af/Vf 

where M/N I means that the terms in the product are arranged from left to right in the 
order of decreasing values of M/N while M/N |" implies opposite ordering of the factors, 
and Uy are defined as in (4.3) with Q(7;t a ) replaced by Q^. 

Noting that the operators for different k > 1 commute, one may combine them into 



00 



Vy = H U k y = exp iJ2 e*y > «(T) = m ~ 2 n(l/m) , (4.6) 

k=l \e=l / m|7 

and rewrite (4.4) into a product over primitive vectors only, 

A 71 , 72 = J] Vy. (4.7) 

7=M7i+iV72, 
M>0,7V>0,gcd(M,iV)=l 



Using the definition of Q given in §2.2 the wall-crossing formula takes the form 

II ^M7i+V 7 2 = H ^Af7i+V 7 2 ' (4-8) 

M>0,N>0, M>0,JV>0, 
gcd(M,N)=l,M/Nl gcd(M,AT)=l,M/JVt 

where are defined as in (4.6) with Q replaced by 0, ± . 

The invariants Cl~(M7i + ./V72) on one side of the wall can be determined in terms of 
the invariants VL + {M^fi + iV~7 2 ) on the other side by expressing both sides of (4.8) into single 
exponentials using the Baker-Cambell-Hausdorff (BCH) formula, and then comparing the 
coefficients of each e 7 on either side. These equations can be solved iteratively to determine 
Cl~(M , -fi + iV7 2 ) in terms of the f2 + 's, starting with the lowest values of (M,N). This is 
most conveniently done by projecting the relation (4.8) to the finite-dimensional algebra 

A m ,n = A/{ R 

e m 7 i+n72 } , (4-9) 

m>M and/or n>N 
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and using the Baker-Campbell-Hausdorff (BCH) formula to commute the factors (see later). 
Since Am,n is a finite dimensional algebra generated by 77171 + 7172 for < m < M and 
< n < N, we have a finite number of equations relating (77771 + 7772) to Q + (77771 + 7772). 
For example a trivial consequence of (4.8) is the relation 

IT (Af 7 i + iV 72 ) = fi + (Af 7 i + iy 72 ) (4.10) 

whenever M = or N = 0. This follows from the fact the algebras Am,o and Ad.tv ar e 
commutative. 

In order to derive the semi-primitive wall-crossing formula and generalizations thereof, 
it is also practical to work with the infinite dimensional algebra Am,oo, an d consider the 
generating functions 

00 00 
Z ± (M, q) = J2 fi± ( M 7i + N l2) q N , Z ± (M, q) = ^ ^(M^ + N l2 ) q N , (4.11) 

N=0 N=0 

for fixed value of M. These two objects are related by 

Z±{M,q) = ^^(M/d,/) , Z ± (M,q) = J2^-Z ± (M/d,q d ) . (4.12) 

d\M d\M 

4.2 Charge conservation from KS formula 

In this section we shall draw attention to one specific feature of the wall-crossing formula 
given in (4.8), namely 'charge conservation'. It follows from the algebra (4.1), and the 
definition of given in (4.6) that after combining each side into a single exponential, the 
coefficient of e 7 consists of a sum of products of the form rL^(7«) with X^7« = 7' U P to 
an overall numerical constant. Thus, any relation that follows from (4.8) has the property 
that the sum of the charges in the argument of Si's has the same value for all the terms on 
either side of (4.8). Thus when we solve this to find ^"(7) in terms of ^l + { r y) and products of 
f2 + (7i)'s, each term in the expression will have the property that the charges in the argument 
of Cl + in the product will add up to 7. This is precisely the 'charge conservation' rule that 
followed from the use of Maxwell-Boltzmann statistics in §2.4. Note that the wall-crossing 
formula written in terms of fi ± does not have any such manifest charge conservation. We 
shall see examples of this 'charge conservation' rule in the explicit examples to be described 
below. 

When several walls are crossed consecutively, the black hole molecules and bound molec- 
ular clusters can be decomposed into smaller molecules and eventually just single atoms. This 
is the attractor flow conjecure [6]. Knowing the indices of the atoms, one can in principle 
determine the contribution to the index of the total molecule, or equivalently flow tree. 
These flow trees are naturally parametrized by nested lists, e.g. ((71, 72), 73), which need 
to satisfy the 'charge conservation' rule. The structure of the nested lists is identical to the 
commutation relations of the KS formula in terms of ^(7), and allows to determine easily 
the contribution of a molecule to the index [28]. 
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4.3 Primitive wall-crossing 

In A\ t i the BCH formula reduces to 

e X e Y = e X + Y + \[X,Y ] ; (413) 

since all multiple commutators involving three or more generators vanish. The wall-crossing 
equation takes the form 

exp(fi + (7i)e 7l ) exp(fi + (7i + 7 2 )e 7l+72 ) exp(fi + (7 2 )e 72 ) 

= exp(0-(7 2 )e 72 )exp(0-(7! + 7 2 )e 7l+72 )exp(0"(7 1 )e 7l ) . (4.14) 

From this we find the primitive wall-crossing relation 

n- (i, i) = n+(i, i) + (-i)^ 7l2 n+(o, i)n+(i, o) , (4.15) 

where for simplicity we have denoted f2 ± (M7i + A^7 2 ) by f2 ± (M, N) and (71, 72) by 712. 

4.4 Generic 3-body and 4-body contributions 

We shall now extract the generic n-body contribution from the KS wall-crossing formula. To 
explain what is meant by 'generic', let n = 3, and a±, a 2 , 03 be three distinct (not necessarily 
primitive) elements of T such that their central charges Z ai , i = 1,2,3 are arranged in 
clockwise order in c_. Then = (a*, otj) > for % < j. For definiteness we shall choose the 
ctj's such that a i2 > a 23 . In this case, in the convention described below (3.1), the different 
linear combinations of the a^s will follow clockwise order as we move from left to right in the 
list (3.8). Furthermore, as we move from the left to the right in this list, the central charges 
will follow clockwise order in the chamber c~ and anti-clockwise order in the chamber c + . 
We can now ask the following question: what is the coefficient of Q + (ai)Cl + (a 2 )Q + (a 3 ) in 
the expression of fl~ (a.\ + a 2 + a 3 ) — Q + (ai + a 2 + a 3 ) in terms of sum of products of Cl + 's ? 
We refer to this coefficient as the generic 3-body contribution to wall-crossing. In order that 
the KS formula be consistent with the explicit computation of bound state degeneracies of 
black hole molecules, this coefficient must agree with the quantity g(a\, a 2 , 03) computed in 
§3. 

In order to carry out this computation, we can pretend that all Q + (a) other than those 
for a = a\,a2,oi3 vanish. In this case the left hand side of (4.8) takes a simple form 
V+ 3 V+ 2 V^. Our task is to manipulate this product so that in the final expression V^'s follow 
the ordering of the en's given in (3.8). For this it is most expedient to expand each in 
(4.8) in a Taylor series expansion in ± (a), manipulate the left hand side so that in each 
term the e Q 's follow the order (3.8), and then identify the coefficients of each term on two 
sides. Since our goal is to find the term involving f2 + (ai)f2 + (a 2 )f2 + (a3), we can focus on 
the term Q + (ai)fl + (a 2 )fl + {a 3 )e a3 e a2 e ai in the Taylor series expansion on the left hand side. 
Repeated use of (4.1) gives 

e as e a2 e ai Co 1 Ca 2 Ca3 ~\~ ^((0^2; < -^l))6o. 1 +a 2 CQ3 ~\~ K ( {oi 3 , Ct\ ) )6 Ql _|_ Q3 6 Q2 

(4.1b) 

+ «((a 3) a 2 ))e ai e a2+a3 + re((a 2 , ai))«((a 3 , "1 + o; 2 ))e ai+a2+a3 . 
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Note that the e^'s on the right hand side follow the order given in (3.8). Since in the Taylor 
series expansion on the right hand side of (4.8) e ai+a2 +a 3 is multiplied by fl~(a\ + a 2 + #3), 
by comparing the coefficient of e ai+a2+a3 in the left and right hand side of (4.8) we get, from 
(4.16), 

AQ(ai + a 2 + a 3 ) = (-l^+^+ais a 12 (« 13 + a 23 ) fi + (a 2 ) Q + (a 3 ) + ■■■ (4.17) 

where • • • represent terms other than the one containing the product Cl + (ai) Cl + (a 2 ) ^ + («3)- 
Happily, this agrees with the result (3.11) from the black hole bound state analysis. 

In a similar fashion, let us consider the generic 4-body contribution. We assume the 
same ordering of the different linear combinations of the a^s as given in (3.13). Then one 
finds for the jump across the wall 

AO(«i + a 2 + a 3 + a 4 ) = - + («i) + (a 2 ) + (a 3 ) + (a 4 ) x 

[k(«i 2 )k (a 13 + a 23 ) k (atu + a 2i + "34) 

+k(oi 13 )k (aiu + a 34 ) k (a 2 i + a 23 + a 24) (4-18) 
+K(a 23 )K(a 14 )K ((a 2 + a 3) «i + a 4 ))] + . . . 
= (-l) 1+E ^ Q ^ + (cn)fi + (a 2 )0 + (a 3 )fi + (a 4 ) x 

[ai2 «13 «24 + «13 «14 «24 + «12 «23 «24 + «14 «23 «24 

+a i2 a 23 a 34 + a 13 a 23 a 34 + a 14 a 23 a 34 + a 13 a 24 a 34 ] + . . . , 

where the dots represent other contributions. This again agrees with the result (3.15) of 
the black hole bound state analysis. For completeness also give the result for the two body 
contribution to Af2: 

Ati( ai + a 2 ) = (-l) Ql2+1 ai 2 Q + { ai ) ti + (a 2 ) + • ■ ■ . (4.19) 

Finally note that this method can be easily generalized to the case when some of the 
Qii's are equal. For example if we are looking for a term Yii (^ + ( a i)) m ' i n the expression 
for Cl~ (^2 i m i ai), then we must expand the left hand side of (4.8) so that in V ai we keep 
the {pt + {oii) e ai ) m * /mil term, then carry out the rearrangement described above, and finally 
identify the coefficient of e^\ m ^ Qi in the resulting expression. Alternatively we could simply 
use the generic n-body formula for non-identical particles, take the limit when several of the 
ccj's approach each other and then include the symmetry factor 1/ JT^ mj in accordance with 
(2.16). It is easy to see how this rule arises from the KS formula, - given a factor of {ep) m for 
any vector (3 we can replace it by eg(i)e«(2) • • • e«( m ) for m distinct vectors /3^\ . . . , (3^ with 
slightly different phases, 15 carry out the abovementioned manipulations for rearranging the 
vectors and then take the limit when all the /3^'s approach each other to recover the desired 
result. There is a priori an ambiguity in this procedure since in the final configuration the 
relative ordering between two vectors which only differ by the replacement of (3^ by (3^ 
for some pair is arbitrary, but this does not affect the final result since changing this 
relative order picks up a commutator factor that vanishes as (3^ — > (3^\ Thus the only 
effect of having identical particles will be the Boltzmann symmetry factor l/m\. 

15 For this manipulation we can ignore the fact that the charges are allowed to take values on a lattice. 
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4.5 Semi-primitive wall-crossing formulae and generalizations 

A general wall-crossing formula involves computing f2~ (77771 + 7172) in terms of Q + (k / -fi +£72) 
for k < m, i < n. We define the order of the wall-crossing formula as the smaller of m and 
n. In this and the following subsections we give wall-crossing formula? for increasing order, 
starting with order one in this section and ending at order three. For simplicity we shall give 
the result for the case when m is fixed to be 1, 2 or 3, but the result can be easily generalized 
to the case when n is 1, 2 or 3 (see §4.6). In the D6-D0 example described in appendix B, 
order corresponds to the number of D6-branes or the rank of the sheaf. 

4.5.1 Order one 

To extract the semi-primitive wall-crossing formula from the KS formula, we project (4.7), 
(4.8) to the algebra Ai tOQ : Thus we have 

^71,72 = ^yi ^71+72 ^71+272 ^72 = ^72 " ' ^71+272 ^71+72 ^71 • (4.20) 
Noting that e 7l+ 7v 72 all commute in -4.i j00 , this can be rewritten as 

A 7li72 = e x i e Y = e Y e x i (4.21) 

where 

00 00 00 

xf = ^(71 + n 12 ) e 11+Nl2 , y = Y. 0+ (^ 2 ) e^ 2 = E °~(^) e *» ■ ( 4 - 22 ) 

Af=0 t=\ 1=1 

It follows from (4.21) that 

Xf = e- y X+e y . (4.23) 
To evaluate the right hand side of (4.23) we first observe that, for a single term in Y, 

00 _ 

e -o +( *72)e, 72 g7i+7V72 e n+(,72)e, 72 = £ [ K (t ll2 ) fi^)]* e ll+(N+lni)l2 • (4.24) 

Thus, 



ni=\j 



1 ™l 



Y ^/tt K^12)0 + (^72)]' , 

e e 7l+A r 72 e = ^ M J_ — ] e 7l+(A r +E< lni)l2 , (4.25) 

{n k } \ t l ' 



and from eq, (4.23), 



1 "c 



AT = t V(-n + N, 2 ) £ n ^f^' ) e 71+(iv+E , fc(h5 . ,4.26) 

N=0 {n k } \ t e ' J 

Using the relation between Xf and f2 _ (7i + X72) given in (4.22) we now get 

N 

n _ (7i + ^72) = Yl + N 'l*) fi haio(7i, iV - N') (4.27) 

7V'=0 
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where we defined 



^haio(7i,A0= Yi n 



{n k } i 
J2kn k =N 



(4.28) 



In terms of the partition function Z ± (l, q) = Z (1, q) defined in (4.11), we obtain 



Z (l,q) = Z + (l,q) Z halo ( 7l ,g) 



(4.29) 



where 



^haio(7i, q) = E JWTi, = exp ^(-1)<^ 2 > ( 7l , £ 72 > fi+(^ 72 ) 



(4.30) 



N=0 



reproducing (2.22), (2.23). 



4.5.2 Order two 



We now extend the semi-primitive wall-crossing formula to order 2, i.e. compute Z (2,q). 
To this aim we project (4.7), (4.8) to the algebra ^.2,00, 



.4 



71:72 



nV + V + 

^7i+fc72 ^2 7 i- 



■(2fc+l) 72 



.fe=0 



• v + = v~ ■ 

72 72 



v: 



27i+(2fc+l)72 71+^72 



.fc=oo 



(4.31) 



We can combine all factors of V^ +N and V 2 ~ 1+N72 on either side into a single exponential, 
by using the level two truncation (4.13) of the BCH formula. Thus, we now have 



^71,72 e 



(4.32) 



where X±, Y are the same as in (4.22), while 



+N72 



N=0 



Qf(N) = n ± (2 ll + N l2 )±j <\l -^712) ^(Ti + iTa) " ± (7i + JI2) • 



i+j=N 



Eq.(4.32) now implies that 



(4.33) 



(4.34) 



XI = e~ Y Xfe Y X- = e~ x X^e Y . 



(4.35) 



Using the analog of (4.27) with 7l — > 2 7l and eq, (4.33) we get 



A 1 



^(N) = Vj(N') fi halo (2 7l , N - N') 



N'=0 



(4.36) 



- 38 - 



where f2haio is defined in (4.28). Combining this with (4.27), (4.34), we arrive at 
fT(2 7l + iV 72 ) = E ^ + (27i + iV72) ^10(271,^-^) 



0<N'<N 



+ \ E «((j-07i 2 )fi + (7i + ^2)fi + (7i+J72)fihaio(2 7 i,iV-z-j) 

0<i<j;i+j<N 

~2 E -«)7l2)fi + (7l +«72)^ + (7l +j'72) ^halo(7l^ - 0^halo(7l)J ~ j') 

0<j<i,i+j=N 
0<i'<i,0<j'<j 

(4.37) 

This result generalizes eq. (4.10) in [43]. 

The partition functions (4.11) for M = 2 are most conveniently expressed in terms of 
"modified partition functions" 



as 



N=0 

Z ± (2,q) = Zf{q) ± g) * g) 



(4.38) 
(4.39) 



where the star product is defined by 



E art* * E M = E - J')7i2) at hi q i+ i . (4.40) 
i=o i=o o<j<i 

The relation (4.36) then simplifies to a simple wall-crossing identity for the modified partition 
functions, 

Z 2 -(g) = Z 2 + (g)Z halo (2 7l ,g) (4.41) 

where Z halo (27!,g) is given by the same formula as in (4.30) with 7l — > 2 7l . Thus, the 
effective description is still in terms of a halo of non-interacting Boltzmannian particles 
around a core with effective degeneracy Qf(N). 

4.5.3 Order three 

The projection of (4.7), (4.8) to the algebra v4.3 !00 reads 



.4 



71,72 



nv + v + v + v + 

V li+k~i2 ^37i + (3fc+l) 72 ^2 7 i+(2fc+l) 72 ^3 7l + (3fc+2) 72 



.k=0 



■ K 



~;2 



v. 



12 



II ^3 7 i + (3fc+2) 7 2^2 7l + (2fc+l) 72 ^3 7l + (3fc+l)72 ^71+^72 



,fc=oo 



(4.42) 



We can combine factors according to the 7 i-charge by using the level three truncation of the 
BCH formula: 

N 1 1 N 

l g(e*°e^...e*") = E^+9 E [X i ,X j ] + -^[X^X,]] 



i=0 
1 

+ 6 



0<i<j<N 



i,j=0 



£ + 

\0<i<j<k<N 0<k<j<i 



E )[X h [X j ,X k }} + .... (4.43) 

<j<i<N/ 
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Thus, we now have 

A 7li72 = e x i+ x i +x t e Y = e Y e^+^+^a" . (4.44) 
Xf,Xf,Y are the same as in (4.33), (4.22), (4.21) while 



Xf = ^3 (37i + N l2 ) e 3yi+N , 2 (4.45) 

N=0 

where we defined 

Qf(N) =£^(37! + #72) 

± \ E 2i l 7i2) ^(71 + ^2) ^(271+^72) 

+ 7^ E - 07i2)] 2 [fi ± (7i + n 2 )] 2 ^ ± (7i+J7 2 ) < 446 ) 

2j+j=A r 

+ ^ K (*' fc ) n± ^ + ^ 2 ) ^ ^ + ^ 2 ) n± ^ + ^) , 

0<«<j<fc 



«(i,j,A;) = - j) ll2 ) + j - 2k) 112 ) + -K((k - j) 112 ) + k - 2i) 112 ) . (4.47) 


Note that (4.47) may be written differently using the Jacobi identity 

-k(0'-z')7i2) K{{i+j-2k) 112 )-K{{k-j) 112 ) K((j+k-2i) ll2 ) + K ((fc-i) 7l2 ) K((i+fc-2j) 7 i 2 ) = , 

(4.48) 

and that j, k) is symmetric under exchange i -H- fc. As before, eq, (4.44), besides leading 
to eqs.(4.23) and (4.35), gives 

X 3 - = e~ Y Xte Y . (4.49) 
Using the analog of (4.27) with 71 — > 371 and eqs.(4.45) we get 

5^(JV)= E ^(^0^0(371,^-^) (4.50) 

0<N'<N 
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where fihaio is defined in (4.28). Putting together (4.46) and (4.50), we arrive at 
ft- (3-ft + N l2 ) = ^ + (37i + N' l2 ) fi halo (3 7 i, N - N') 

N> 

+ \ E ^(|j-2^|7i2)|ft + (7i+^2)ft + (27i+J72)fihaio(3 7 i,iV-^-j) 

*>0,j>0 
i+j<N 

+ 7^ E ^ ~ *)^ 2 )] 2 ^(jl + + ^ fi halo(3 7 l, iV - 2i - j) 

i>0,j>0 
2i+j<N 

+ E «(^J,A;)ft + (7i + z7 2 )ft + (7i+j7 2 )ft + (7 1 + A;72) fthalo(37i, N — i — j — k) , 

0<i<j<fc 

+ ^ E 2Z ' 712 ) 0+ ^ 71 + ^ 0+ ^ 71 + J ' 72 ) fi halo(7l, i ~ 0fthalo(27lJ ~ /) 

i+j=N 
0<i'<i,0<j'<j 

+ \ E 2 ^ 12 ) K d Z - ^l7i2)n+( 7 i + i'72) ^ + (7i + ^2) n + ( 7 i + / 72 ) 

0<i'<i;fc+Z<J 

fthaio(7i> » - «')fthaio(27i, j -k-l) 

2i+j=iV 
0<i',j"<j,0<j'<j 

fthaio(7i> i - ^haio(7i ; i - i") fthaio(7i> j - f) 
E ^ fc ) + ^2) ^ + (7i + J'72) ^ + (7i + k' l2 ) 

i+j+k=N;i>j>k>0 
0<i'<i;0<j'<j;0<k'<k 

fthalo(7l) « - fthalo(7l>j - /) ^halo(7l, & ~ &') 

+ ^ E K (b' + ^ - 2i|7ia) «(| j - ^|7i 2 )ft + (7i + i'72) ^ + (7i + fc'72) ^ + (7i + i'72) 
o<i'<'t,o<j'<i,o<fc'<fc 

fthalo(7l) « - ^halo(7l> J - j') ^halo(7l) fc ~ fc ')- 

(4.51) 

As before, we may introduce the "modified partition functions" 

00 

Zi(q) = Y,ni(N) q N , (4.52) 

7V=0 

in terms of which the wall-crossing relation takes the simple form 

Zz(q) = Z+(q)Z halo (3 ll ,q) , (4.53) 

where Zhaio(37i, q) is defined in (4.30). The modified partition functions are related to the 
original partition functions Z^(3, q) by a relation of the form 

Zf{q) = Z ± {3, q) + Z ± (l, q) * ± Z ± (2, q) + Z ± (l, q) * ± Z ± (l, q) * ± Z±(l, q) (4.54) 

with the definitions of the 'star' products following from (4.46). 
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4.6 712 > case 

We shall now briefly discuss what happens when 7 12 > 0. We can of course use all the 
formula? derived in this section with 71 -H- 72, but then e.g. eq, (4.51) will give the wall- 
crossing formula in the charge sector N^i + 372. So if we want to find the wall-crossing 
formula for 371 + ^72 we cannot get the result by exchanging 71 and 72 in (4.51). Instead 
we keep 712 > and carefully examine how the subsequent equations are affected. It is easy 
to see that the only place where the sign of 712 enters is in eq, (2.6); for 7 12 > the < sign in 
(2.6) is replaced by a > sign. So if we continue to define Q as in §2.2 then Q + will denote 
the index in the chamber in which the multi-centered bound states of black hole molecules 
exist. We would however like to define Q + as the index associated with single black hole 
molecules, and for this we exchange the definitions of Q + and Q~ . Thus for example in eq, 
(4.51) we have to now exchange Q + with Q~ so that we have an expression for Q + in terms 
of Q~ . We can in principle solve these equations iteratively to find fl~ in terms of Q + , but 
we shall now suggest a simpler method. For this note that exchanging Q + and Q~ in the 
wall-crossing formula (4.8) is equivalent to changing the order of the products on both sides 
of (4.8). This in turn is equivalent to keeping the same order as in (4.8) but changing the 
sign of all the structure constants in the algebra (4.1). This can be achieved by changing 7 12 
to —712- Thus if we replace 712 by — I712I in the formula? we have derived (e.g. (4.51)), we 
shall get the correct wall-crossing formula for both signs of 712. 

4.7 Refined wall-crossing and motivic invariants 

We have already introduced the refined invariants in §2.6. In this subsection we shall review 
the KS motivic wall-crossing formula that computes the jump in the refined index across walls 
of marginal stability and compare it with our wall-crossing formula based on the analysis of 
supergravity bound states. 

In order to describe the motivic generalization of the wall-crossing formula [10], we 
consider a set of generators {e 7 ,7 G T} satisfying the quantum torus relations 

e 7 ey = (— yY 1,1 ' e 7+ y . (4.55) 

The associated Lie algebra is 

K M = (i-y) M) - (-y)- { ™' } ) e 7+y . (4.56) 

Let us also introduce the quantum dilogarithm, 

{xy) n 
„ )...(l 

n=0 v a ' v " ' n=0 

This satisfies the pentagon identity 

E(xi)E(x 3 ) = E(x2)E(x 12 )E(x!) (4.58) 

where 07, £2 ar e two non-commutative variables satisfying x\Xijy = yx2%\ = — X12, and 
reduces to the ordinary dilogarithm for logy = v — y 0, 

B(x, = exp (-i- L, 2 (*) + ^— )V + S + ...). (4.59) 
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We attach to the charge vector 7 the generator 16 

% = ]^ E (y n e 7 r ( ~ 1)nfW7) . (4.60) 

The motivic version of the KS wall-crossing formula [10, 33, 34] again states that the product 

^71,72= II ( 4 - 61 ) 

7=^71+^72, 
M>0,N>0 

ordered so that as we move from left to right the corresponding Z 7 's are ordered clockwise, 
stays constant across the hyperplane of marginal stability ^(71,72)- 

As in the classical case, it is advantageous to combine the generators £4 7 for k > 1 
in a single factor V 1 . For this purpose, rewrite the operator f/ 7 , using the identity for the 
quantum dilogarithm 

EM-xpg^M^), (4.62) 

as 

f/ 7 = exp \ \ k ( y k_ y -kf n j ■ (4.63) 

Then the product of t/g 7 over all £ > 1, 7 being a primitive vector, can be written in terms 
of f2 re f(7,?/) defined in (2.38): 

% = H = ex P ( £ ^ref (^7, !/) ev 7 ) , e 7 = _ 1 . (4.64) 
This is the motivic generalization of (4.6). The wall-crossing formula now takes the form 

II ^M7i+^72 = II Vmji+Nj 2 ' (4.65) 

Af>0,JV>0>0, M>0.V>0>0, 
gcd(M,N)=l,M/Nl gcd(Af,JV)=l,M/JVt 

where V ± are computed using the fi re f(7, y) in the chambers c . It follows from (4.56), 
(4.66), (4.64) and (4.65) that expressed in terms of Q re f( , y,y) the wall-crossing formula will 
satisfy manifest 'charge conservation laws'. 

To see how this reduces to the classical KS formula (4.8) in the limit y — > 1, note that 
the generators 

e 7 = lim e 7 , e 7 = - — - (4.66) 

y -*i y — y 

satisfy the commutation relations (4.1). Moreover in this limit Q Te f( , y,y) approaches ^(7). 
Thus V defined in (4.64) reduce to V ± and we recover (4.8). 



16 In supersymmetric gauge theories, for a vector multiplet with unit degeneracy, U 1 reduces to J7 7 



E (ye 7 ) E (y 1 e 7 ) while for a hypermultiplct one has C/ 7 = E (e 7 ) 



-2 
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4.8 Semi-primitive refined wall-crossings and its generalizations 

The rescaled generators e 7 = e 7 /(y — y^ 1 ) satisfy the same Lie algebra as (4.1), 

[e 7 ,ey] = K((^,-f'),y) e 7+y , (4.67) 
provided k(x) is replaced by its quantum deformation 

K(x,y) = ^~ V ^ ~^J^ = (-l) x sinh(^x)/sinh^, i/ = lnj/. (4.68) 

Moreover, the generators V 1 in (4.64) can be obtained from their classical counterpart (4.6) 
by replacing £l(y) (->■ f2 re f(7, y), e 7 \-> e 7 . Therefore, the wall-crossing formulae derived in 
Section 4 and Appendix A carry over to the motivic case by just replacing 

fi(7) H> fi rcf (7,y) , Q(7) h-> fi ref (7,y) , k(x) h-)> . (4.69) 

In particular, the primitive wall-crossing formula takes the form [15, 33] 

AO ref ( 7 71 + 72, y) = { ~ y> ' ~ , — Orf (7i, V) ^cf (72, 2/) (4.70) 

y- l/y 

while the refined semi-primitive wall-crossing formula is given by 

Z~(l,q,y) = Z + (l,q,y) Z hSLlo (y u q, y) (4.71) 

where 

^haio(7i, 9, 2/) = exp ( — _ -1 ^ref(^72, y) <f ) ■ (4.72) 



=1 



On the right hand side Sl re f (£72, 2/) can be computed in either chamber. (4.72) is in perfect 
agreement with (2.42), showing that the Boltzmann gas picture correctly reproduces the 
semi-primitive motivic wall-crossing formula. To compare (4.72) with known results, we 
note that in terms of the "integer" motivic invariants Q TC f, n {l), (4.72) becomes an infinite 
product [33] 

00 fc|7ia| 

z^ (jx,q,y) = n n n ( x - (-i)*-^^*- 1 -^)™*^™ . (4.73) 

k=l j=l neZ 

To see this, note that the logarithm of the r.h.s. of (4.73) can be rewritten as 

"EEE £'^(-1)^0^72) [(-l)^- h ^y n+J »- 1 - k ^\] d . (4.74) 

d>l k>t n j=l 

The sum over n leads to ^lr e {,n{k^f2,y d ), while the sum over j is geometric, leading to 

bg z halo ( 7 i, q, y ) = E -a n «*(*^» y d ) d _ V — ? M > ( 4 - 75 ) 

d>i k>i y y 
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where we have used 712 < to replace I712I by —712- Setting i = kd, the sum runs over 
divisors d of N and reproduces (4.72). 

The order two and three motivic wall-crossing formula can be obtained by making the 
replacements (4.69) in eqs.(4.37), (4.51). 

Finally, let us consider the problem of determining a generic 3-body contribution to the 
wall-crossing formula: given three charge vectors cti, a 2 and a 3 in T, what is the coefficient 
of f2+ f (ai,y)Q^ e{ (a 2 , 2/)^ref( a 3' V) 1X1 t ne expression for Af2 re f (ai + a 2 + a 3 , y)l The analysis 
is a straightforward generalization of that in §4.4 and the final result is obtained by replacing 
K{otij) by K,(aij,y) in (4.16). This gives: 



This is in perfect agreement with (3.10) computed from the spectrum of bound states of a 
3-centered configuration in supergravity. Similarly the 4-body contribution can be computed 
by replacing k(o) by n(a,y) in (4.18). The result is in perfect agreement with (3.14) and 
(3.51). We have also carried out a similar computation for 5-body contribution and compared 
with the results obtained by following the procedure of §3, but we shall not give the details. 

4.9 KS vs. supergravity 

Eventually one would like to prove that the KS wall crossing formula given in (4.65) is 
equivalent to the one obtained from quantization of multi-black hole solutions as given in 
(1.5), (1.6). We have not yet reached this goal, but would like to point out some common 
aspects of these two formula?. The summand in (1.6) depends analytically on the ay's, 
but the analyticity of the sum is broken by the third condition described below (1.6). In 
particular this constraint measures whether ^2 a P^ a \ represented as a vector in the two 
dimensional plane in the convention described below (3.1), lies to the left or the right of the 
vector a.\ + ■ ■ ■ + a n . Let us denote by B the set of all vectors of the form £\ m^ai with m, = 
or 1, and arrange them in an order following the convention described below (3.1). Let B' 
denote the subset of elements of B which lie to the left of the central element ai + • • • + a n . 
In this case the expression for <7 re f given in (1.6) depends on the subset B' . As we vary the 
«j's this subset may change and in that case g rc { will be given by a different analytic function 
of the a^-'s. Note however that g rci does not depend on the relative ordering of the vectors 
inside the subset B'. 

Now this lack of analyticity is also manifest in the KS wall-crossing formula. To see this, 
recall the procedure for manipulating the KS formula given in §4.4. Here we are supposed 
to begin with the product e an ■ ■ ■ e Ql and bring it to the order in which the vectors appear in 
the set B. Changing this order leads to a different final order of the e^'s and hence we expect 
5 l ref(tti, • • • , ttn, u) i given by the coefficient of e ai+ ... +oln in the final state, to change. This in 
turn prevents g re f to be given by an analytic formula involving the ay's for all choices of a^. 

Now, the KS prescription for computing g TC f a priori seems to depend on more informa- 
tion than (1.6) since the KS formula requires the detailed ordering of the vectors in B, rather 
than just the unordered list of the ones which lie to the left of ai + • • • + a n . We shall now 



AQ(ai + a 2 + a 3 ,y) 



(_ 1 )«i2+a 23 +«i3 (ginhi/)- 2 S inh(i/ai 2 ) smh(u(a 13 + a 23 )) 
ti + (a 1 ,y)ti + (a 2 ,y)ti + (a 3 ,y) + -- - . 



(4.76) 
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show that the KS formula in fact only depends on the unordered list of vectors which lie to 

the left of ai + h a n . For this let us consider a given order of all the a^'s and supppose 

that we have brought e an ■ ■ ■ e ai to the required order. Now consider the effect of switching 
the relative order between two neighbouring vectors (3i and fa on the left of ai + • — h a n . 
This will require to switch the corresponding ep 1 and ep 2 and will produce an extra factor 
of e/3 1+/ 3 2 . But since fa and fa both lie to the left of a\ + ■ ■ ■ + a n , fa + fa will also lie to 

the left of «i H + a n . Thus such switchings can never produce a factor of e ai+ ... +OCn . The 

same argument holds if we switch two vectors on the right of a,\ + • • • + a n . Thus a term 

proportional to e QlH hQ , n can arise only if we switch a vector from the left of a,\ + . . . a n with 

a vector to the right of a\ + ■ ■ • + a n . This shows that the non-analyticity of (1.6) and the 
KS formula are controlled by the same data. 

5. Wall-crossing from the Joyce-Song formula 

In their work on Donaldson-Thomas invariants for coherent sheaves on a Calabi-Yau three- 
fold X [14] (which presumably count D6-D4-D2-D0 bound states in type IIA string theory 
compactified on X), Joyce and Song give a fully explicit expression for the rational DT 
invariants Q~ on one side of the wall, in terms of the rational DT invariants Q + on the other 
side. Thus, the JS formula can be viewed as the solution to the implicit relation given by KS. 
In particular, it directly provides the functions #({«;}) appearing in (2.16), i.e. the solution 
to the black hole bound state problem. It should be noted however that the JS wall-crossing 
formula involves sums over many terms with large denominators and large cancellations, 
and is less computationally efficient that the KS formula (compare for instance table (5.16) 
with the analogous computation for KS given in eq. (4.18)). In addition, the simple rule for 
dealing with identical particles mentioned at the end of §4.4 is not at all obvious from the 
JS formula. 

One way of understanding the large redundancy is that JS work with Abelian categories, 
where constituents are either a subobject or a quotient object of the complete object. In 
physical terms this means that different terms in the JS wall-crossing formula keep track of 
the order in which the constituents (molecules) make the complete object (bound molecular 
cluster). But in physics (and in the derived category on which KS analysis is based) such 
a distinction is not present. For example the existence and index of a bound state of two 
primitive constituents carrying charges 71 and 72 is independent of the order in which we 
choose 71 and 72. For this reason, the JS formula contains many terms which must combine at 
the end to ensure the independence of the final answer of the order in which the constituents 
are chosen. The KS formula is less redundant, at the cost of being implicit and perhaps less 
rigorously established. 

After stating the JS formula in §5.1 and its implication for the black hole bound state 
problem in §5.2, we verify the equivalence of the JS and KS formulae for generic three and 
four body contribution, and semi-primitive wall-crossing in §5.3 and §5.4. 

5.1 Statement of the JS formula 

In [13, 14], the authors define rational- valued generalized Donaldson-Thomas invariants 
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DT (r) for any class 7 G C(X), where X is a Calabi-Yau three-fold, C(X) is a positive 
cone 17 inside K{X) and r is a stability condition. They furthermore establish a general 
wall-crossing formula for the variation of DT 7 (r) under a change of stability condition. 
Conjecturally, the rational invariants DT 7 (r) are related to integer invariants DT 7 (r) by 
a relation identical to (1.2), and DT 7 (r), DT 7 (r) are particular instance of the invariants 
0(7; t), $7(7; t) considered in [10]. The stability condition r is determined by a point in 
Kahler moduli space i. In the convention we have chosen, 



In this section we shall assume that the conjectured relation between the rational and integer 
invariants holds, and denote the rational DT invariants of [13, 14] as 0(7; i). The JS wall- 
crossing formula then furnishes the solution to the KS formula (4.5), i.e. expresses ^"(7) in 
terms of 

Let ai, a 2 , • • • , a n be n charge vectors in the positive cone C{X) inside the charge lattice 
described by eq, (2.5) , and (t,t) be a pair of points on the Kahler moduli space with 
associated stability conditions (r, f). To express the JS wall-crossing formula, we first need 
to introduce two functions S(ati, . . . , a n \ t, t) and U («i, . . . , a n \ t, t), whose role is to capture 
the relevant information about the ordering of phases of {Z a .}. We define S(ai, . . . , a n \ t, t) G 
{0, ±1} as follows. If n — 1, set S(ai, t,t) — 1. If n > 1 and, for every i — 1, . . . , n — 1, 
either 

(a) r(ttj) < r(a i+ i) and f{a\ H h oti) > f(a i+1 H h a n ), or 

(b) r(aj) > r(a i+ i) and f(«i H h otj) < r(« i+ i H V a n ) , (5.2) 

let S(ai, . . . , a n ; t, t) = (— l) r , where r is the number of times option (a) is realized; otherwise, 
S(a u . . . ,a n ;t,t) = 0. 

To define U(cxi, . . . , a n ; t, t), consider all ordered partitions of the n vectors a« into 
1 < m < n packets {a nj l+ i, • • • , a aj }, j = 1, . . . , m, with = a < ai < • • • < a m = n, such 
that all vectors in each packet have the same phase r(aj). Let (3j = a aj _ 1+ i + • • • + a aj , j = 
l,...,m be the sum of the charge vectors in each packet. We refer to the ordered set 
{(3j,j = 1, . . . ,m} as a contraction of {ai}. Next, consider all ordered partitions of the m 
vectors f3j into 1 < I < m packets {/3& fc _ 1+ i, ■ • • , /\J, with = bo < b\ < ■ ■ ■ < bi = m, 
k = 1, . . . , /, such that the total charge vectors 5k = /3& fc _ 1 +i + ■ • • + /3& fe , k = 1, . . . , I in each 
packets all have the same phase f(5k) (which is then equal to ^(X] a i))- Now associate to 
each of the I packets in the contraction a factor S(f3b k _ 1+ i, /3& fc _ 1+ 2, . . . , t, t) as defined 
above, and define the ^/-factor as the sum 



r{l,t) 



argZ(7,t) . 



(5.1) 



(-1) 



rn 



1 



U(oii, a n ; t, t) = 



,...,(3 bk ;t,t) •] 



(aj - aj_i)! 



(5.3) 



over all partitions of ctj and partitions of /3j satisfying the conditions above. 



17 This is the analog of the wedge f introduced in eq. (2.5). 
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Finally, departing slightly from the presentation in [13], let us define the £ factor 



£(«!,... ,a n ) = ^ (o!i»«i>» ( 5 - 4 ) 

connected oriented trees: edges i — 

vertices { 1 , . . . ,n} , edge i — yj, implies i < j 

where the sum runs over all connected trees g with n vertices labelled from % — 1 to % = n. 
We denote by = {1, . . . ,n} the set of vertices, and by g^ 1 ' the set of oriented edges (i, j), 
with the orientation inherited from the standard order i < j on g(°\ In other words given 
any labelled tree, and an edge of this tree connecting i to j, we associate to this edge an 
orientation from i — y j if i < j. In order to implement this formula on a computer, it is 
useful to note that there are rS n ~ 2 ^ labeled trees with n vertices, which are in one-to-one 
correspondence with their Priifer code, an element in {1, . . . ,n} n_2 . 18 

Having defined the S, U and £ factors, we can now state the JS wall-crossing formula 
(eq. (5.13) in [13]): 

- ~ ^ x ^ (-l^-i+E^IH,^)! „ " _ 

^(li t )=}_^ 2^ o^=I U(a 1 ,...,a n ;t,t)£(a 1 ,...,a n ) \[n(ai,t) , 



n>l ( ai ,...,a„)£C(X) i=l 
7=«iH hOn 



(5.5) 



where the second sum runs over all ordered decompositions 7 = ot\ + • • • + a n with ai G 
C(X). Note that due to this constraint, eq, (5.5) is automatically consistent with the charge 
conservation property. 

The JS formula is valid for any pair of points in the moduli space with stability conditions 
r, f. Now we restrict to the vicinity of a wall of marginal stability "P(7i, 72)- The only states 
whose BPS invariants can jump are those whose charges lie in T C C(X). Moreover, their 
discontinuities only depend on the BPS invariants of states in Y. We denote the states 
7 = Mji + iV72 by 7 = (M, N). As in previous section, let us assume that 712 < 0, and take 
t in the chamber c + where multi-centered configurations are absent, and t in the chamber 
c~ where they are present. From (2.6), we see that, in the vicinity of the wall, 

±(7, 7 ') [r(7,c ± )-r(7 / ; C± )] >0 . (5.6) 

Therefore, the conditions (5.2) become 

(a) (cci, a i+ i) < and (a± H \-ai,a i+ i-\ h a n ) < 0, or 

{b) a i+ i) > and (ati + • • • + a iy a i+ \ + • • • + a n ) > , (5.7) 

and (5.5) gives an expression for Q~ in terms of Cl + . 

18 Note that these labeled trees differ qualitatively from the attractor flow trees. One can view these labeled 
trees as connecting the cndpoints of the attractor flow trees. 
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3- 




Figure 1: The 16 labelled trees contributing to C(cti, a 2 , a^, a^). 



To illustrate the afore going prescription, we evaluate in detail the term proportional to 
^ + (7i) 2 ^ + (72) 2 in the expressions for f2 _ (27! + 272). We first evaluate the U factors, 



^(72,72,71,71) 


= 1^(272,271) + ^ 


5'(72,72,27i) + -5(272, 71, 71) + 5(72,72,71,71) 


^(72,71,72,71) 


= 5'(72,7i,72,7i) - 


^(72, 71) 2 = 


^(72,71,71,72) 


= -^(71,72)^(72 


7i) + ^(72,271,72) + 5(72,71,71,72) = 


^(71,72,72,71) 


= -^(71,72)^(72 


7i ) + ^5(71,272,71) + 5(71,72,72,71) = 


^(71,72,71,72) 


= 5(71, 72, 71,72) - 


^5(7i, 72 ) 2 = - 


^(71,71,72,72) 


= ls(2 7l ,2 72 ) + l 


5(7i,7i,272) + ^5(271,72,72) + 5(71,71,72,72) 



For example in the first line the only non- vanishing contribution comes from the 5(27 2 , 2ji)/4 
term. For the other terms the condition (5.7) fails for at least one i. To compute the C 
factors we observe that there are 16 trees with 4 labelled nodes (see Figure 1). Out of those, 
twelve are obtained by various inequivalent permutations of the tree connecting nodes 1 to 
2, 2 to 3 and 3 to 4, and four are obtained by inequivalent permutations of a tree that 
connects node 1 to each of the nodes 2, 3 and 4. The ^-factors are computed by adding the 
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contributions from each of these 16 trees, leading to 

£(«!, Oi 2 , Ol 3 , Qt 4 ) ="l2"l3"l4 +«12 023«14 + "l3"23"l4 + ai3"24"l4 

+ "23 "24 "l4 + "l2 "34 "l4 + "23 "34 "l4 + "24 "34 "l4 
+ "12 "13 "24 + "12 "23 "24 + "13 "23 "24 + "12 "13 "34 
+ "l2 "23 "34 + "13 "23 "34 + "12 "24 "34 + "l3 "24 "34 



(5.9) 



This gives 



£(72, 72, 7i, 7i) = -47i 3 2 , £(7a> 7i> 72, 7i) = 27 3 2 

£(72,7i,7i,72) = , £(71,72,72,71) = (5.10) 
£(71,72,71,72) = -27? 2 , £(7i,7i,72,72) = 47i2 

The total contribution from these terms to the right hand side of (5.5) is thus given by 

n-(27i + 2 72 ) = \ll 2 t^ + (7i)] 2 [^ + (7 2 )] 2 + • • • , (5.11) 
in agreement with the formulae (4.37) and (A. 2). 

5.2 Index of supersymmetric bound states from the JS formula 

It is useful to rewrite the JS formula (5.5) as a sum over unordered decompositions of 
the charge vector 7, 

7=«lH ha n 

where |Aut({aj})| is defined as follows. We choose any specific order of the aj's and identify 
|Aut({aj})| as the order of the subgroup of the permutation group S n on n elements which 
leaves invariant this particular ordering of the {a^j's. In particular if all the aij's are different 
then the automorphism subgroup is trivial and |Aut({ttj})| = 1. On the other hand if m of 
the aj's are the same then |Aut({aj})| — ml. The coefficient g({ai}) is then given by 

{_l\n- 1 +'Ei <j (°>i, a j) 

g{{oii}) = 2^ U (ov(i), a a ( n y,t+,t-) C (a a{1 ), . . . , ov( n )) , 

(5.13) 

where the sum runs over all permutations of n elements. Relabelling the set {a{\ into 
a partition {m rjS } such that 7 = (M,N) = ^2 rs {r, s)m r s , the coefficient g{{a.i}) is the 
same as the one appearing in (2.16). Thus it should be identified with the index of the 
supersymmetric quantum mechanics of n distinguishable charged dyons in R 1,3 , along the 
lines of [7]. 

5.3 Generic 2-body, 3-body and 4-body contributions 

Let us first derive the primitive wall-crossing formula 7 — > 71 + 72 from the JS wall crossing 
formula. In this case there is only one tree - connecting the nodes 1 and 2 with the arrow 
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directed from 1 to 2. If we choose the first node to be 71 and the second node to be 72 then 
since 712 < possibility (a) in (5.7) is realized, and we have {7(71,72) = S (71,72) = — 1- 
Since n = 2, (5.13) now gives a contribution of ^(— l) 1+7l2 7i 2 (— 1). An identical contribution 
comes from the term where we put 72 in the first node and 71 in the second node since now 
we have U = S = 1 and the (a a ^,a a ^) factor in (5.13) is now equal to 721. Adding the 
two contributions we recover the primitive wall-crossing formula (4.15). We summarize the 
S, U and £ factors for the two permutations in the table below. 



(7(12) 


S 


U 


c 


12 


a 


-1 


712 


21 


b 


1 


-712 



(5.14) 



Next we reproduce the result of Sec. 4.4 for the generic three body contribution to the 
wall-crossing from centers carrying charges (a\, a 2 , olz) . The order of the a^s is given by 
(3.8). With the definition of the ordering explained below (3.1) we see that if (a,/3) follow 
a clockwise order then (a, (3) > 0. Since the phases are assumed to be generic, the U and 
S factors coincide, i.e. U(ai, atj, a^.; t, t) = S(ai,aj,ak]t,t), for any permutation {i,j,k} of 
{1, 2, 3}. In this case, there are three trees contributing. The S, {7, £-factors are summarized 
in the table below. Substitution of these data into eq. (5.5) reproduces directly Eqs. (3.11), 
(4.17). 



(j(123) 


S 


u 


£ 


123 


bb 


1 


ai2«13 + «13«23 + «12«23 


132 


b- 





«12«13 - «13"23 - «12«23 


213 


ab 


-1 


-«12«23 + «13«23 - ai2«13 


231 


-a 





ttl2«13 - «13«23 - «12«23 


312 


ab 


-1 


ttl3tt23 - ttl2«23 - ttl3«12 


321 


dbQb 


1 


«13«23 + «12«13 + Oii 2 a 2 3 



(5.15) 



Next, we turn to the generic 4-body case. We use the same phase ordering as eq. (3.13). 



The non-vanishing contributions are given in the following table 



cr(lZo4 ) 


c 
o 


U 


C 


1 OQ A 

1264: 


bbb 


1 


"l2"l3"l4 +ai2tt23«14 + "l3 "23 "l4 


+ "13 "24 "14 








+"23 "24 "14 + "12 "34 "14 + "23 "34 "14 


+ "24 "34 "14 








+"l2 "13 "24 + "12 "23 "24 + "13 "23 "24 


+ a±2 CK13 «34 








+«12 «23 "34 + "13 "23 "34 + "12 "24 "34 


+ "13 "24 "34 


1 O A O 


bba 


— 1 


"l2 "13 "14 - "l2 "23 "14 ~ "l3 "23 "l4 


— "13 "24 "14 








+"23 "24 "l4 + "l2 "34 "l4 ~ "23 "34 "l4 


— "24 "34 "14 








— «12 "13 "24 + "12 "23 "24 + "13 "23 "24 


+ "12 "13 "34 








-"12 "23 "34 - "13 "23 "34 ~ "12 "24 "34 


- "13 "24 "34 


1 /I OQ 


bab 


— 1 


«12 «13 «14 + «12 «23 "14 + "13 "23 "14 


— "13 "24 "14 








— a 2 3 "24 "l4 — "l2 "34 "l4 ~ "23 "34 "l4 


+ "24 "34 "14 








-Q!i2 "l3 "24 - "l2 "23 "24 ~ «13 "23 "24 


- a 12 a 13 a 34 








— "12 "23 "34 — "13 "23 "34 + "12 "24 "34 


+ "13 "24 "34 


1 /I oo 


baa 


1 


"12 "13 "14 — "12 "23 "14 ~ "13 "23 "14 


— "13 "24 "14 








+"23 "24 "14 — "12 "34 "l4 + "23 "34 "l4 


+ "24 "34 "14 








— "12 "13 "24 + "12 "23 "24 + "l3 "23 "24 


- "12 "13 "34 








+a 12 a 2 3 «34 + "13 "23 "34 + "12 "24 "34 


+ "13 "24 "34 




bba 


— 1 


— "12 "13 "14 + "12 "23 "14 + "13 "23 "14 + "13 "24 "14 








—"23 "24 "14 + "12 "34 "14 ~ "23 "34 "14 


— Q; 2 4 "34 "14 








+"12 "13 "24 — "12 "23 "24 ~ "l3 "23 "24 + "l2 "l3 "34 








— "12 "23 "34 — "13 "23 "34 — "12 "24 "34 


— "13 "24 "34 


O /I 1 o 


bab 


— 1 


"12 "13 "14 + "12 "23 "14 ~ "13 "23 "14 


— ai3 Q;24 CKi4 








— «23 Q;24 0(14 ~ "12 "34 "14 + "23 "34 "14 


+ Q;24 "34 "14 








— "12 "13 "24 — "12 "23 "24 + "13 "23 "24 


+ "12 "13 "34 








+"12 "23 "34 - "13 "23 "34 + "12 "24 "34 


- "13 "24 "34 


O /I Q 1 


baa 


1 


— "12 "13 "14 + "12 "23 "14 + "l3 "23 "l4 


+ "13 "24 "14 








— "23 "24 "14 — "12 "34 "14 + "23 "34 "14 


+ "24 "34 "14 








+"12 "13 "24 — "12 "23 "24 ~ "13 "23 "24 


- "12 "13 "34 








+"12 "23 "34 + "13 "23 "34 + "12 "24 "34 


+ "13 "24 "34 


O 1 O /I 


abb 


— 1 


— "12 "13 "14 — "12 "23 "14 + "l3 "23 "l4 


— CK13 «24 "14 








— «23 "24 "14 + "12 "34 "l4 — "23 "34 "l4 


+ "24 "34 "14 








—"12 "13 "24 — "12 "23 "24 + "13 "23 "24 


- "12 "13 "34 








— "12 "23 "34 + "13 "23 "34 + "12 "24 "34 


— "13 "24 "34 


Q1 /lO 

ol4z 


aba 


1 


— "12 "13 "14 —"12 "23 "14 + "l3 "23 "l4 


+ "13 "24 "14 








+"23 "24 "14 + "12 "34 "l4 — "23 "34 "l4 


— "24 "34 "14 








+"12 "13 "24 + "12 "23 "24 ~ "l3 "23 "24 


- a 12 a 13 a M 








-"12 "23 "34 + "13 "23 "34 ~ "12 "24 "34 


+ "13 "24 "34 


3241 


aba 


1 


— "12 "13 "14 —"12 "23 "14 ~ "13 "23 "14 


+ "13 "24 "14 








+"23 "24 "14 + "12 "34 "l4 + "23 "34 "l4 


— "24 "34 "14 








+"12 "13 "24 + "12 "23 "24 + "l3 "23 "24 


+ a 12 a 13 a u 








+"12 "23 "34 + "13 "23 "34 ~ "12 "24 "34 


— "13 "24 "34 
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ct(1234) 


s 






4213 




1 


"12 "13 "14 ~t~ "12 "23 "14 — 


"13 "23 "14 ~\~ "13 "24 "14 








"1-CK23 "24 "14 — "12 "34 "14 + 


"23 "34 "l4 — "24 "34 "l4 








T"12 "13 0X24 T "12 "23 "24 


"13 "23 "24 "12 "13 "34 








-\-OL\1 «23 «34 "l3 "23 "34 


s~\> s~\> s~\> J s~\ , s~\> s~\> 

"12 "24 "34 ~r "13 "24 "34 


4321 


fl fl fl 


_1 


— «12 «13 tti4 — «12 "23 "14 _ 


"13 "23 "14 - "13 "24 "14 








— «23 "24 "14 — «12 "34 "14 ~ 


"23 "34 "14 — "24 "34 "14 








— "l2 "l3 "24 — "12 "23 "24 ~ 


"l3"23«24 - "12 "13 "34 








-«12 a 2 3 "34 - "13 "23 «34 ~ 


Oli2 OL24 "34 — "13 "24 "34 



Adding up all these contributions, one finds again the result of Eqs. (3.15), (4.18). 
5.4 Semi-primitive wall-crossing formula from JS 

Let us now derive the semi-primitive wall-crossing formula from (5.5), i.e. compute ^(7; c_) 
in terms of ^(7; c + ) for 7 = (1, N) G T. We mostly follow the discussion in [41], Section 3, 
suitably generalized. 

At given order n, the most general ordered decomposition 7 = J2i=i n a * * s 



7 = (0, N x ) + ■ ■ ■ + (0, jVi,_i) + (1, A^) + (0, N im+ i) + ■ • • + (0, N n ) 



(5.17) 



where {Ni} is a partition of N of length n, Yli<i<n^i = an( ^ ^ — ^* — ^ l a bels the 
position of the charge (1, 0) in this decomposition. For such a pair ({Ni}, one may check 
that the 5" factor is given by 



^({AVKi*; c+,c_ 



-1^ 



(5.1? 



if = 1 or 2* = 2, and vanishes otherwise. To compute the U factor, note that due to 
the phase constraints, the ordered partition of {ai} labelled by {a,} must decompose into 
ordered partitions of the 'head' set {Ni}, i — 1, . . . , — 1 and 'tail' set {Ni}, i — i* + 1, . . . , n, 
while the vector (l,Ni t ) must lie in its own packet. Moreover, since the phases of all the 
vectors /3j, except the vector (l,N it ), are the same in the chamber c_, all the /3j's must be 
be grouped into a single packet at the second step in order to satisfy the phase constraint on 
the SkS. Otherwise the packet containing (1, A 7 ^) will have different phase from the others. 
The ?7-factor therefore reduces to 



U({Ni},i*;c + ,C- 



i=i 



l<m<n 
0=ao<ai<...a m =n 



a,- - a 7 _i ! 



(5.19) 



Now, let be the packet in which the vector (l,AjJ lies. In view of (5.18), the factor 
S((3i, (3 2 , ■ ■ ■ , /3 m ; c + , c_) vanishes unless = 1 (which happens if z* = 1) or = 2 (which 
happens whenever > 1 and we group all the vectors (0, Ni), . . . , (0, A^-i) in a single 
packet). In either case, the contraction of the head set is trivial, and we are left with 



U({Ni},i*;c + ,C- 



—y 

- 1)! ^ 



_iyen(p) 



(5.20) 
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where the sum runs over ordered partitions of the n — 2* elements in the tail set, i.e. integer 
sequences {p r , 1 < r < len(p)},p r > 1 such that ^2 r p r = n — i*. This evaluates, for all 
> 1, to the binomial coefficient 

U( { N i} ,i,;c + ,c.)= {h (5.21) 

a result which is in particular independent of the partition {N}. 

Now, we turn to the sum over graphs g. Due to the Landau factors (a iy aj) in (5.5), the 
only contributing graph is a tree rooted at z*, with leaves 1, . . . , z* — 1, z* + 1, . . . , n. The 
Landau factor is then given by 

i* — 1 n 

£= n<M) 5 (i,^j> n (mj^m)). (5.22) 

Taking into account the additional factor l/2 n_1 , we thus arrive at 



n— i* _,n— 1 



v ' > 2 n ^ 1 (z* - 1)\ n 



n>l Ki»<n 



(5.23) 



Plugging this expression in the partition function (4.11), we can easily carry out the sum 
over Ni and obtain 

n>l Ki»<n v * 7 v * 7 



where Zh a io(7i, <?) is the same function introduced in (4.30). The sum over z* leads to 

halo (7l ; < 

(n-1)! 



z- M = z+(i, q) £ N^fri.dr 1 (525) 



n>l 

and the sum over n finally leads to 

Z~(l,q) = Z + (l,q)Z halo (li,q) ■ (5.26) 

Finally we note that the derivation given above can be simplified using the Boltzmann 
gas picture in which we analyze identical particle contribution to wall-crossing as a limit 
of non-identical particle contribution and then include an extra symmetry factor 1/N\ for 
N identical particles. To see how this works, we consider the case where we have (N + 1) 
different c^'s satisfying > for i < j, with the understanding that we shall eventually 
take the limit in which the first N aij's approach 72 or its multiple and a/v+i approaches 71. 
Now the S and the U factors coincide as in §5.3. Furthermore since eventually we shall take 
the limit in which the first N a^s coincide, the requirement of a non- vanishing C tells us 
that only trees which contribute are those in which a/v+i is connected to all the other cc^'s. 
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We can still choose the direction of the arrows arbitrarily. Let us consider a configuration 
in which m of the arrows are directed towards a^+i and n = N — m are directed away from 
it. In this case the arrows go from the nodes 1 to m towards the central node (m + 1) and 
from the central node to the nodes (m + 2) to (N + 1). The Landau factor associated with 
these nodes is (— l) m YliLi( a N+i, ou)- Next we need to assign the aj's to the different nodes. 
First of all there are ways of deciding which of the a^s will be assigned to the first m 
nodes. Once this is done there is no further freedom of rearranging the aj's among the first 
m nodes or the last (N — m) nodes; in order that U = S does not vanish, the aj's must 
be arranged in a clockwise order among the first m nodes and anti-clockwise order among 
the last (N — m) nodes. In this arrangement the possibility (a) is realized (N — m) times 
and the possibility (b) is realized m times. Hence the corresponding U is given by (—l) N ~ m . 
Substituting these into eq, (5.13) we now get 



N /N\ N N 
g({a t }) = 2- JV (-l)^<>^> E Uj II K+i,^) = (-1)^<>^> Y[(a N+1 , ai ) . 

m=0 ^ ' t=l i=l 

(5.27) 

If we now take the set (ai, . . . , ajv) to contain mi copies of lj2 with ^Imi = N then, 
including the symmetry factor Yli l/ m ^5 we S e t the coefficient of the f2 + (7i) Yli^ + (h2) mi 
term to be 

II^tK- 1 )' 712 ^)]^, (5.28) 



which is the desired result. 



Acknowledgement: We are grateful to D. Joyce and J. Stoppa for useful correspon- 
dence. J.M. thanks the LPTHE for hospitality during part of this work. The work of J.M. 
is partially supported by ANR grant BLAN06-3-137168. A.S. acknowledges the support of 
Chaires Internationales de Recherche Blaise Pascal during his stay at LPTHE where the 
work began, of the J. C. Bose fellowship of the Department of Science and Technology, India 
and of the project 11-R& D-HRI-5.02-0304. 



A. Wall crossing formulae in special cases 

In this appendix we give explicit wall-crossing formula? in some special cases. These cases 
illustrate the general results of §4 and §5, and serve as tests of the equivalence of the JS and 
KS wall-crossing formula? . For brevity we state the results for the classical rational invariants 
f2, but the same formulae also hold for the motivic rational invariants with the replacement 
(4.69). For notational convenience we shall denote fi ± (M7 1 + A r 7 2 ) by fi ± (M, N). The result 
for (M,N) = (2,2) below is in agreement with [28], Eq (2.13). 

Afi(l,2) =/<27i 2 )fi + (0,2) fi + (l,0) 

+ \n ll2 ]f [fi+(o, i)] 2 n + (i, o) + nfa) n + (o, i) n + (i, i) (A ' 1} 
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Afi(2, 2) = «(47i 2 ) H + (0, 2) fi+(2, 0) + «(2 7l2 ) 0)fi + (l, 2) + fi + (0, l)fi + (2, 1)] 

+ «( 7 i2) «(27i 2 ) ^ + (o, i)n + (i, o)n + (i, i) + \nii2)} 2 «(2 7l2 ) + (o, i) 2 n+(i, o) 2 (A . 2) 

+ ^(2 7l2 )] 2 p + (2, 0)fi+(0, l) 2 + n+(0, 2)fi + (l, 0) 2 ] 



Afi(l, 3) = «(3 7 i 2 ) H + (0, 3) 0) + «(2 7 i 2 ) H + (0, 2) 1) + K ( 7l2 ) Q + (0, 1) 2) 

+ K ( 7l2 ) k (2 7i2 ) n + (o, i) n + (o, 2) o) + i[ K ( 7l2 )] 2 n + (o, i) 2 i) 



+ ^K 7 i 2 )fn + (o I i) 3 n + (i ) o) 



(A.3) 



Aft(2, 3) = «( 7l2 ) i) 2) + «(3 7 i 2 ) o) 3) + «(6 7i2 )q + (o, 3) n+(2, 0) 

+ «;(4 7l2 )ft + (0, 2) fi + (2, 1) + /<2 7l2 )n+(0, 1) fi + (2, 2) 

+ ^H7i 2 )] 2 fi + (o, 1) n+(i, 1) 2 + ^H3 7l2 )] 2 n + (o, 3)n + (i, o) 2 

+ \ [K712)] 2 + [k(2 7i2 )] 2 + «( 7l2 ) «(3 7l2 )] fi+ (0, 1) 0) 2) 

+ ^ [«( 7l2 ) k(2 7i2 ) + «( 7l2 ) k(4 7i2 ) + k(2 7i2 ) k(3 7i2 )] n + (0, 2)ft+(l, 0) 1) 

+ K (2 7l2 ) «(4 7l2 ) fi + (0, 1) ft + (0, 2) ft + (2, 0) + ^(2 7l2 )] 2 fi + (0, l) 2 n+(2, 1) 

+ \ (3[«:(7i 2 )] 3 + «(37i 2 ) [^(7i 2 )] 2 + «( 7 i2) K2 7l2 )] 2 ) [O+(0, l)] 2 0) n+(l, 1) 

+ ^( 7l2 ) k(2 7 i 2 ) [k( 7 i 2 ) + k(3 7 i 2 )] 0)] 2 fi + (0, 1) fi+(0, 2) 

+ ^(2 7l2 )] 3 [fi+(0, l)] 3 n+(2, 0) + ^K 7l2 )] 3 [3«( 7 i 2 ) + «(37i 2 )] ^ + (0, l) 3 n + (l, 0) 2 

(A.4) 
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Afi(2,4) 



1 

12 



1 

48' 



k(27i 2 )k(7i2) 4 + 7^(4 7 i 2 )k( 7 i 2 ) 4 



n + (o,i) 4 n+(i,o)< 



+— [4/ t (2 7l2 ) K ( 7l2 ) 3 + K (4 7l2 ) K ( 7l2 ) 3 + /€(2 7l2 ) 3 K ( 7l2 )] n+(o, i) 3 n+(i, o)n+(i, i) 



+ 



-/t(2 7l2 ) 2 fi:(7 12 ) 2 + -/€(2 7l2 )K(4 7l2 )fi:(7 12 ) ; 



ft+(0,l) 2 fi+(0,2)ft+(l,0) ; 



1 



+ i fi + (0,l) 2 n + (l,l) 2 K (2 7l2 )/s:(7 12 ) 2 



+ 



- K (27 12 ) K ( 7l2 ) 2 + -/€(4 7l2 ) K (7 12 ) 2 + -«(2 7 



12 J 



n + (o,i) 2 n + (i,o)n + (i,2) 



+ [k(2 7i2 ) 2 /< 7i2 ) + k(27 12 )«(4^ 12 )«( 7i2 )] fi+(0, 1)0+ (0, 2)0+(l, 0)O+(l, 1) 
+ K (2 7l2 ) K ( 7l2 )n + (0, l)fi+(l, 2) + ^(4 7l2 ) 2 fi+(0, 4)fl + (l, 0) 2 

+ l - [k{2 1i2 )k{ 1i2 ) + «(4 7i2 )k( Ti2 ) + «(2 7l2 )«(3 Tl2 )] ft + (0, l)fi+(l, 0)O+(l, 3) 

+ l - [k(2 7i2 )/<3 7i2 )k( 7i2 ) + k(3 7i2 )k(4 7i2 )/< 7i2 )] ft+(0, l)ft + (0, 3)0+(l, 0) 2 

+ ^ [«(67 12 )k( 7i2 ) + k(2 7i2 )k(37 12 ) + «(3 7i2 )k(4 7i2 )] fi+(0, 3)D + (1, 0)O+(l, 1) 

+ ^-K(2 7l2 ) 4 n + (0, l) 4 fi + (2, 0) + ^(2 7l2 ) 3 O + (0, l) 3 fi+(2, 1) 
24 o 

+ ^(2 7 i 2 ) 2 n + (0, 2)fi+(l, l) 2 + ^(2 7 i 2 ) 2 n + (0, l) 2 fi+(2, 2) 

+^(4 7l2 ) 2 n+(o, 2) 2 n+(2, o) + «(2 7l2 )o+(i, i)n + (i, 3) + /<2 7l2 )o + (o, i)n+(2, 3) 

+ ^(27 12 ) 2 K (47 12 )n + (0, 2) 2 fi+(l, 0) 2 + ^K(27 12 ) 2 /s:(47 12 )n + (0, l) 2 ft + (0, 2)fl+(2, 0) 

+k(4 7i2 )0+(1, 0)O + (l, 4) + /<4 7l2 )fi+(0, 2)0 + (2, 2) + /<6 7l2 )O + (0, 3)Q + (2, 1) 
+ K (2 7l2 ) K (4 7l2 )n + (0, 2)fi+(l, 0)fl + (l, 2) + «(2 7l2 )«(4 7l2 )ft + (0, l)ft + (0, 2)fl+(2, 1) 
+k(2 7i2 )/€(6 7i2 )O + (0, l)n + (0, 3)0+(2, 0) + «(8 7l2 )fi+(0, 4)0+(2, 0) . 

(A.5) 



To derive these results from the JS wall-crossing formula described in §5, we use results 
for the U and £ factors tabulated below. In these tables, the first column describes the total 
charge, the second column lists the allowed ordered decompositions and the third and the 
fourth columns give the U and £ factors introduced in §5.1. Note that the U factor typically 
involves a sum of multiple S-factors representing possible partitioning of the constituents as 
in (5.8), while the £ factor comes from a sum of several trees as in (5.9). 
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B. D6-D0 bound states 



In this subsection, we test and apply the wall-crossing formula on generalized DT-invariants 
for dimension zero sheaves on a Calabi-Yau threefold X, for which many results are already 
known in the literature [10, 42, 41, 51]. 

The stability conditions for coherent sheaves on X depend on the complexified Kahler 
moduli t a = B a + iJ a . The holomorphic central charge Z 7 is given, in the large J limit, by 
Z-y = — J x e~ t ^Td(X) A 7. Let us denote by f2(r, n; t) the generalized DT-invariant for a 
sheaf of rank r, vanishing first and second Chern class, and third Chern class n. Physically, 
Q(r,n) counts the number of bound states of r D6-branes and n DO-branes, with charge 
vector 7 = r7i + n7 2 , 712 = — 1. 

It is known that for infinite volume and small £>-field, a configuration of r > 1 D6-branes 
and n > 1 DO branes do not form any bound state [52]. Moreover, there are no bound states 
of r > 1 D6-branes, while n > 1 DO bind into precisely \x\ bosonic (fermionic) bound states 
for negative (positive) x> where x is the Euler number of X . Thus the only non- vanishing 
DT invariants in this chamber are [10, 13] 

fi + (l, 0) = 1 , fi+(0, n) = - X (n > 0) , (B.l) 

where we have used the notation n ± (m,n) = (77171 + 71,72). Similarly, the motivic invari- 
ants are given by 

n+f(l,0,y) = 1 , ^ r + c f(0,n,i/) = -P(y)/y 3 (n > 0) , (B.2) 

where P(y) = 1 + &2£/ 2 — b^y 3 + 6 2 |/ 4 + y 6 is the Poincare polynomial of X, such that (B.2) 
reduces to (B.l) in the classical limit y — > 1. 

By increasing the magnitude of the 5-field, one reaches the wall of marginal stability 
^(71,72). We refer to the chamber across this wall as the 'DT' chamber. We shall obtain 
the motivic DT invariants in this chamber for r < 3 by applying the formulae derived in 
Section 4, suitably generalized to the motivic case according to the discussion in Section 4.7. 

To apply the semi-primitive wall-crossing formula? and its higher order generalizations, 
it is useful to introduce partition functions 

00 00 
Z±(M, q,y) = J2 n tef(Mli + N l2 , y) q N , Z±(M, q, y) = £ n± (M 7 i + iV 7 2, y) q N , 

(B.3) 



N=0 N=0 



for fixed value of M. These two objects are related by 



Z^(M, q,y)=Y: "f* Z±(M/d, q d , y d ) 

d\M Vy y ' 



d\M yy y > 



y d ) 



In the chamber c + , the modified partition functions (B.3) are simply 

2(1+7) ' z - f (3,g,2/) = - 3(] 



(B.4) 



y 2 



Z+ (1, q, y) = 1 , Z+ (2, g, y) = — — , Z+ (3, q, y) = » ■ ^ . (B.5) 
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Moreover, the partition function of the halo degeneracies follows from (4.73), 

oo k 

z^ {i^y) = fin{(i- (-qfy^T 1 t 1 - (-riV"*)" 6 " (i - (-^V^f 3 

k=l j=l 

(1 - (-q) k y k ~ 2j r b2 (l-i-qfy'-^y 1 } 

(B.6) 

where we have made a change of variable j — > k + 1 — j in the last two terms to make the 
y —7- y~ l symmetry manifest. Applying (4.71), we immediately obtain the partition function 
of the motivic invariants with r = 1, 

Z~ { (1, q, y) = (1, q, y) = Z ha io(7i, q, y) . (B.7) 

In the classical limit y — > 1, this reduces to 

Z-(l, ? ) = Z-(l, ? ) = [M(-g)]*, (B.8) 

where M(q) = YYk=iO- ~ q k )~ k is the Mac-Mahon function [53, 54, 55, 56]. For comparison 
with the higher rank formula? of [42, 41] below, it is useful to note that the expansion of this 
formula in powers of \ reads 

oo p p oo oo , 

e -xE£,Ho g (M- 5 « = 1 + I: X nEE|(-9) fen< ' < B£! > 

p=l i=l ki = l rii = l 1 



For r = 2, we have 



U 



Z-<(2, q, y) = Z- f (2, q, y) - » Z~ f (1, g 2 , y 2 ) . (B.10) 
We now use (4.34), (4.41) together with the replacement (4.69) to obtain 



Z rcf (2, q, y) = 2 {i + y 2 [^10(271, g, y) - Zhaio(7i, <? 2 , y 2 )] 



7 X] K d ni _n 2|,l/)fi rcf (l,ni,y)fi rcf (l,ra2,i/)g 



ni+n 2 



(b.ii; 



4 

ni,n 2 



In the classical limit, this reduces to 

Z"(2, g) = V(g)] 2 * - \[M{ -g 2 )] x - ^ ("ir^K - *a| n _ (l,m) fi-(l,n 2 )g' 



ni+ri2 



4 

ni,ri2 



(B.12) 

This agrees with [42],Thm 1.2 and [41], (2. 9), who obtain 

^(2,n)=-5„, + ^2- 2 ^ £ Il^-i^^-/ 2 ) 

p=l ^ fci>0,n,>0 i=l 1 

T,ikm=n (B.13) 
~\ zZ {-^) ni ~^\ni-n 2 \n-{l,n x )n-{l,n 2 ) 

ni>0,n2>0 
n\+n2=n 
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where it is understood that the third term is zero if n is odd. 
For r = 3, we use similarly 



Z; ci (3,q,y) = Z- { (3,q,y) - - £ - Z~ f (1, g 3 ,y 3 ) (B.14) 



and (4.46), (4.53) to obtain 



y 2 

Z rel ( 3 ' ^ V) =^7^ YV~A [ Z halo(37l, Q, 2/) ~ ^halo(7l, <? 3 , ?/ 3 )] 

3(1 + y l + y 4 L J 



ni+n2 



- re(|7i 2 - 2n x |,?/)^ ref (l,ni,?/)n ref (2,n2,y)g 

ni,ri2 

4(X + y 2) E K ( 2 \ n 2-ni\,y)fi; c{ {l,n 1 ,y)Q; c{ (l,n 2 ,y 2 )q 



2\ _,ni+2ri2 



ni,n 2 



^ £ «( ni - na, y) 2 (1, m, y)] 2 fi" f (1, n 2 , y) g 2 " 1 ^ 



12 

ni,n 2 



/t(m, n 2 , n 3 , y) ft ref (1, m, y) ft ref (l, n 2 , y) ft ref (l, n 3 , y) g ni_H 

"■1>"2>"3 



(B.15) 



In the classical limit, this reduces to 



Z-{^ q )= l -[M{- q )Y*- l -[M{- q *)Y 



il+n 2 



" ^ E (-l)" 2 " 2n >2 - 2^| m) 0-(2, n 2 ) <f 1 

ni,n 2 

- \ E l n2 - ml n"(l."0 n-(l,na) <f 1+2 ™ 2 (B lg) 

ni,n 2 

E(m-^) 2 [^(l,m)] 2 ^(l,n 2 )g 2 " 1+ " 2 

- ^ K(n 1 ,n 2 ,n 3 )^(l,n 1 )^(l,n 2 )^(l,n 3 )g ni+n2+n3 . 

rai>n 2 >n3 

This can be rewritten as in [41], eq. (2.12), after correcting the coefficient of the 3-body 
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term from 1/4 into 1/6 in that equation: 



E tl k :~y 1 



»! 14 — ' rii 9 

^ fci>0,rii>0 i=l 



lT(3,n) = ^ n ,o + XX 
p=i 

E (-l) n2_2 ">2-2n 1 |Q-(l,n 1 )fi-(2,n 2 ) 

ni>0,ri2>0 
ni+ri2=n 

-| E I 77,1 _ n 2 \Cl~(l,n 1 )Q~(l,n 2 ) 

ni>0,ri2>0 
2ni+n2=n 

E (ni-n 2 ) 2 [n-(l,n 1 )] 2 n-(l,n 2 ) 

ni>0,ri2>0 
2ni+ri2=n 

~\ E [(rii - n 2 ){n 1 + n 2 - 2n 3 ) + (n 2 - n 3 )(2n 1 - n 2 - n 3 )] 

0<ni<ri2<n3 
ni+n 2 +n 3 =n 

Q-(l,m)fi-(l,n 2 )fi-(l,n 3 ) 



(B.17) 



where it is understood that the third term vanishes if n is not a multiple of 3. 

The classical invariants Q~ (r, n) for low values of (r, n) are summarized in the table 
below: 



r\n 





1 


2 


3 


4 







-x 


-x 


-x 


-x 


1 


1 


-x 


|(x 2 + 5x) 


-|(x 3 + 15x 2 + 20 X ) 


^(X 4 + 30 X 3 + I55x 2 + I26x) 


2 








-x 


-|(x 3 + 15x 2 + 20 X ) 


-^(X 4 + 30 X 3 + I19x 2 + I02x) 


3 











-x 


^(X 4 + 30 X 3 + I55x 2 + I26x) 



(B.18) 



r\n 



-X 

1 - mX (X 4 + 50 X 3 + 575 X 2 + 1630 X + 624) 

2 -y^X (7x 4 + 250 X 3 + 1925 X 2 + 3890 X + 1248) 

3 -i2oX (7x 4 + 250y 3 + 1925 X 2 + 3890 X + 1248) 



(B.19) 



For low values of (r, n), the motivic invariants can also be expressed directly in terms of the 
Poincare polynomial of the CY threefold X, 

(-y) 6 Q- f (l,2,y) J- [P(y 2 ) + [P(y)] 2 ] + (y* + y") P(y) 

= 1 + (6 2 + l)y 2 - 6 3 y 3 + ^ (&i + 56 2 + 2) - (6a + l)^ 5 

+ l -y & (2b\ + Ab 2 + bj-b 3 + 2)- (b 2 + l)b 3 y 7 +\{b\ + 5b 2 + 2) y 8 



b 3 y" + (b 2 + l)y 10 + y 



.12 



(B.20) 
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-y) 9 0- f (l,3,y) =\[P{y)f + \p{y)P{y 2 ) + (y 4 + y 2 ) [P(y)] 2 



+ [^P(y 3 ) + {y 8 + y 6 + y 4 ) P(y) 



(B.21) 



(B.22) 



y ) 13 n; cl (2,4,y) =^[12 (y 3 + y) 2 [P(y)] 3 + {y 2 + 1) [P(y)] 4 

+ 3 (y 2 + 1) (2P (y 4 ) + [P (y 2 )] 2 + 4 (y 8 + y 4 ) P (y 2 )) 
+ 6 (y 2 + 1) (P (y 2 ) + 6y 8 + 4y 6 + 6y 4 ) \P{y)] 2 
+ 4(3 (y 6 + y 2 )P(y 2 )+2((y 2 + l)P(y 3 ) 

+3 (y* + y° + 2y 4 + y 2 + l)y 6 ))P(y) 



(-y) 17 Q- f (3, 5,y) = - (y 4 + l) (6P (y 5 ) + 5P (y 2 ) P (y 3 )) 
+ ^(2/ 4 + 2/ 2 + l) P(y) (2P(y 4 )+P(y 2 ) 2 ) 

+ 1 (y 2 + 1) P(y) ((y 2 + 1) P(y) + 2 (y 6 + y 2 )) p (y 3 ) 

+ ^P(y) P (y 2 ) [(y A + 3y 2 + 1) P{y? + 6 (y s + 2 y « + 2y 4 + y 2 ) p( y ) 

+12 (y* + y 6 + 2y 4 + y 2 + l) y 4 ] 
+ ^P(y) [l20 (y 2 + l) 2 (y 4 + y 2 + l) y 4 P{y) 2 + (y 4 + 5y 2 + l) P{y) 4 

+20 (y 8 + 4y 6 + Ay 4 + y 2 ) P(y) 3 + 240 (y w + 2y 8 + 3y 6 + 3y 4 + 2y 2 + l) y 6 P(y) 
+ 120 (y 12 + y 10 + 2y 8 + 2y 6 + 2y 4 + y 2 + l) y 8 ] 

(B.23) 

The motivic invariants 0~ f (r, n, — y) have an expansion in positive and negative powers 
of y of the form: 

K 
n=-K 

for some integer K with 

c n = c_ n , (-l) n c n > 0, c n e Z , (-l) n (c - c n+2 ) > for n > . (B.25) 

These are required for consistency of the interpretation (2.26) for the refined index. In 
particular (— l) n {c n — c n+ 2) measures the total number of states with angular momentum 
n/2 contributing to the index (without counting the angular momentum factor {n + 1)). 
Thus our results provide strong support to the motivic wall-crossing formula. Furthermore 
in all examples the refined index seems to satisfy the symmetry property 



Q ref (r, n, y) = fi ref (n - r, n, y) 



(B.26) 



when < r < n. We have tested this for all r < 3, n — r < 3. In the classical case, this was 
established in [51]. 
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C. Generalizations of Seiberg-Witten spectra 



Assume that in the region c+, the spectrum consists only of two states 71,72 with BPS 
degeneracies f2 + (7i) = q, ^+(72) = p, respectively. Using the formulae derived in the text, 
we find that the BPS invariants fi~(M7 1 +A^7 2 ) in the region c_ are given, for min(M, N) < 3 
by 



/ jko(-ir ip<i(ip+ 1) \ 

f607 5 p 3 g 2 + 807 4 p 3 ij + 2207 4 p 2 g 2 

— 7P(7P + 1)9 +27 7 3 p 3 + 2407 3 p 2 <j + 2607 3 P9 2 

i(-l) 7 7P9 *7P9(79 + 1) |(-1) T 7 2 P(7P + 1)9 ®yp + M-yp + l)q + 1) +63 7 V + 220 7 2 P9 + IOO7V 

(79 + 1)( 7 9 + 2) (37P + 7 (3 7 p + 2)<j + 1) ( lq + l)( 7 9p + p + 9) (37p + 4 7 (7p + l)q + 2) +42 7 p + 6O79 + 8j 

|7P9 , , - T i 5 7p(7P+l)(7P + 2) 9 

J7P9 (2j(jqp + p + q) i 7 p(7P+l)9 j^"/ 2 p(lP + 1) (7P + 2) (6 + I67P + 87V + 3079 

(79+1) (2 7 p + 3 7 (7P + 1)9 + 1) 9(39 + p(3 7 9 + 2)) +457 2 p 9 + 157 3 p 2 9) 

J(-1) T 7P ij7P(7P+l) i|(5(-l) T 7P(7P+l)(7P + 2) 

9 (~1) T 7P9 J7P<7P+1)9 (7P + !)(7P + 2)9 (7P + 2)(7P + 3)9 (7P + 3) ( 7 P + 4)9 



(c.i; 



where 7 = 712 (In this table, M increases from bottom to top, while N increases from left to 
right). This table has been given for 7 < 0; as discussed in §4.6, for a general sign of 7 the re- 
sult is given by replacing 7 — > — \^y\ in this table. With this replacement the table is symmet- 
ric under p •f-)- q, M •f-)- N, and we can recover some of the results for M = 4 and 5 from the 
corresponding results for N = 4 and 5. For (p, q, 7) = (1, 1, —1) (C.I) agrees with the pentag- 
onal identity (4.58). On the other hand for (p,g,7) = (1, 1, -2), (2, 2, -1), (1, 4, -1), ^+(7) 
gives the strong coupling spectrum of the SU{2) Seiberg-Witten theories with Nf = 0,2,3 
flavors. For example (in a suitable basis) the strong coupling spectrum of the Seiberg-Witten 
theory [3, 4] has a single state of charge 71 = (2, 1) and a single state of charge 72 = (0, — 1) 
for Nf = 0, 19 two states of charge 71 = (1, 1) and two states of charge 72 = (0, —1) for 
Nf = 2 transforming respectively in the (2,1) and (1,2) representations of the SO (4) flavour 
group, and four states of charge 71 = (0, 1) and a single state of charge 72 = (1, —2) for 
Nf = 3 transforming respectively in the 4 and 1 representation of the SO (6) flavour group. 
Thus (C.I) reproduces the weak coupling spectrum of the SU{2) Seiberg-Witten theories 
with Nf = 0, 2, 3 flavors, as demonstrated in [16, 34]. 

If instead the spectrum in c + consists of three states 71,72,71 +72) with respective BPS 



Here (n e , n m ) denotes a state carrying electric charge n e and magnetic charge n. 



degeneracies p, q, r, the spectrum in c_ is given by 



jj79(37(79P + P) 2 
+ 3r .( 7r+ l) 
571(7? + 1) +p(79(2 79 + 3) 

((-l)^p(71+2) + 3r) 3(-l)-'7(579 + 3)r + 1)) 



j7 P9(3(5p97 
+3(p + q)~,+ l)r 
+ 2(-l)7(37r 2 + 3r 
+ (7P+1)(79 + 1) 
(79P + P + 1))) 



f (-1)^79 
f(-l)7p( 7 «j+ l) 
+ 2r) 



(— 1) 7Pi + r 



j7P9 

(2 7 (79P + P 
+4(-l)Tr) - 



J(-1) 7 7P 

(3 7 (57P + 3) 9 r + (-l)7 
Vq (37r 2 + 3r + ( 7 p + 1) 

(-1)7 + 1) 9 (2 7 p + 37(7P+ 1)9 + 1))) 



jj7P(7P + 1) 

I 7 p(7pq + q + 2(-l)7 r ) ((_l)7( 7 p + 2 )q + 3r) 



i7P(9(127 2 (97P + 4)r 5 
+ 12 7 (57P + 2)r 
+ (7P + 1)(3 7 P 
+47(7p+l)<j + l) 
(3 7 P 

+4 7 (7P + 1)9 + 21) 
+ 2f-l)7 r (7(27r 2 ' 
+q(7p(3l7p + 33) 
+7(7p(67 7 p + 93) + 32) 9 
+8)) +4)) 



- 6r 



7 l ?5 7p(_20(-l)7p 2 q(2(7p+ 1) \ 

(7P + 2)q 2 +3r 2 ) 7 4 \ 

+30pr(47r 2 + 6r + ( 7 p + X)q 

(ll 7 p + 7(377p + 32) g + 16)) 

7 2 + 2(-l)7 9 ( 3 07 2 (7p 

(677P + 90) + 30)r 2 

+607(27p + l)(87p + 5)r 

+ (7P + l)(7(207(7p(57p 

(27p+7)+39) + 15)g 2 

+60( 7 p + l)(2 7 p + 1)(27P + 3)l 

+9p(37p(3 7 p + 7) + 14)) + 24)) 

+ 15r(7(7 3 g(1637tj + 75)p 3 

+ 127 2 g(4079 + 13)p 2 

+ (7(4677 3 2 + lllg 

+4r(47r + 9)) + 24)p 

+I5O79 2 + 30q 

+8r(2 7 r + 3)) +8)) 



X57P(7<7 2 9(379 + 2)p 3 
+ 2713(679 + 10(-l) T 7r 
+ 3)P 2 

+ (1579 2 + 4(9(-l)7 7 r + 1)9 
+ 6r(7r + l))p 
+ 16(-l) 7 gr 

+ 6(9 2 +r 2 )) 
-3(-l + (-l)7)r) 



^7p(7p + l)(7p + 2) 
«7P + 3)9 + 4(-l)7 r ) 



J2Q V 1V\ rp ^ 

(57(7p(237p + 55) + 30) 
qr + (-1)7(607(7P + l)r 2 
+ 20(2 7 p+l)r + (7p + 2)9 
( 7 (7P + 2)(8p 
+ 15( 7 p+ 1)9) +6))) 



jij7p(7p + l)(7p + 2) 
(7P + 3)((-l)-'(7P + 4)q 
+ 5r) 



\ p 

(C.2) 

For (p, q, r, 7) = (1, 1, 1, 1) this reproduces again the pentagonal identity (4.58) (read back- 
wards). For (p, q,r,j) = (1,1, 1, —1) it reproduces the weak coupling spectrum of the SU(2) 
Seiberg-Witten theory with Nf = 1 flavor [34], since the strong coupling spectrum of the 
latter has three states of unit degeneracy each, carrying charges 71 = (0,1), 72 = (1,-1) 
and 71 + 72 = (1, 0). 



D. U(N) quiver quantum mechanics from Boltzmann black hole 
halos 

In this appendix we shall show that our Boltzmann gas picture of multi-centered black 
hole system also makes a specific prediction on quiver quantum mechanics without oriented 
loops. Consider a multi-black hole configuration with (n + N) centers, with the first n centers 
carrying charges «i, . . . ,a n and each of the last N centers carrying charge (3q. We take all 
the «j's and (3q to lie in a two dimensional plane, and assume, without any loss of generality, 
that (Xij = (ai,aj) > for i < j. Suppose further that each center carries unit intrinsic 
index, i.e. Q + (cti) = Q + ((3 ) = 1, and that Q + (k(3 ) = for k > 1. In this case the system is 
described by a U (l) n x U(N) quiver quantum mechanics [7], and the refined index Tr (—y) 2j3 
of this quantum mechanics is given by 

(_yyC u(1)nxUiN) ( ai ,...,a n A) P mnxU(N) ( ai , ...,a n , 0o, y) ■ (D.l) 

Here Pu(i) n xU{N){.&i, ■ ■ ■ > Po, y) is the Poincare polynomial associated with a quiver quan- 
tum mechanics with n + 1 nodes, with the first n-nodes carrying U(l) factors and the last 
node carrying a U(N) factor. Furthermore there are arrows from the i'th to the j'th node 
for i < j, \(ai, /?o)| arrows from the ith node to the last node if (aj,/5 ) > and \(ai, f3 )\ 
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arrows from the last node to the ith node if (otj, A) < 0. Cu(i)« x u(N){&i, ■ ■ • , «n, A)) denotes 
the maximum power of in Pf7(i)" X c/(jv)( Q; i) • • • , a n, A; y), signifying the maximum value of 
2 J 3 carried by the system. An explicit expression for the Poincare polynomial of this quiver 
will be given in eq. (D.7) below. On the other hand we can also calculate the index using 
(1.3)-(1.6). First note that we have, from (1.3), 

Using (1.5) we find that this configuration contributes 



E 



<? ref («!, . . . , a n , A x h, 2/3 x fc 2 , . . . , y) H | Lfi+Qfo)* j (D.3) 



to Tr (— y) 2j3 , where g re { is given by (1.6) and the symbol 1/3q x ki means that the vector 
l(3 is repeated ki times. To study the consequences of equating (D.3) to (D.l), note that 
the expression for g rc f given in (1.6), after stripping off the first factor is nothing but the 
Poincare polynomial of a quiver theory with only U(l) factors. This gives 

<7 ref («!,..., a n , A) x fc!,2A xk 2 ,---,y) = (-y)^E, °*-*E<l<o«,A>>l 

X -P(7(l)™x(7(l) fc ixC/(l) fe 2x-( a l5 • • • i a n,Po x fei, 2/3q X fc 2 , • • • , y) ■ 

Using eqs.(D.l), (D.2) and (D.4) we get a direct relation between the Poincare polynomials 
of a quiver carrying a U(N) factor and those carrying only U(l) factors: 

Pu(i)*xu(N){(Xi, ■ ■ - ,(Xn,Po,y) = Y P u(i) n xY\iU(i) k i( a ii ■ ■ ■ > a n, A x h, 2/3 x k 2 ,- ■ ■ ,y) 



1 / n, n -1 \ fci 



(D.5) 



We could also derive a similar relation where the first U(l) n factor is replaced by a product 
of U(M) factors by choosing among the set (a±, . . . , a n ) identical charges. Let us use a 
shorthand notation where Pni,...,n p (Pi, ■ ■ ■ , Pp, y) denotes the Poincare polynomial of a quiver 
with the z-th node having a U(Ni) factor and charge A an d furthermore arrange the A' s 
such that > for i > j. Similarly CN lt ...,N p (Pij • • • > Pp) will denote the maximum power 
of y 2 in this polynomial. Then the generalization of (D.5) will take the form 

y- CN ^"^'-"' M PN 1 ...N p (Pi,...,P P ,y) 

= ^2 PNi,...,N r -i,lxki,lxk2,-,Nr+i,.-,N p (Pl, ■ ■ ■ , A-l> A x 2A x k 2 , ■ ■ ■ ,/3 r+ i,...,P p ,y) 



y 



tt— y ~ y f-r 

(D.6) 
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where we have substituted y — > —y, using the fact that the Poincare polynomial of a quiver 
is an even function of y. By repeated use of (D.6) we can express the Poincare polynomial of 
a quiver with U(N) factors in terms of the Poincare polynomial of quivers with U(l) factors 
only. 

We now test e.q. (D.6) against the Poincare polynomial of a quiver quantum mechanics 
with U(N) factors but no oriented loops computed in [39] (see also [7]): 



P Nl ,..., Np (Pu ■ ■ ■ ,&,y) =(v 2 ~ l) 1 -** y~^ N ^ 

£ (-irv^<-^^ N ^]jm,y]\)^ 



partitions a,i 



(D.7) 



i<j 



where 



v 2N - 1 

i N > y \ - V^T' [N ' y]l ~ [1,y][2, y] "' [iV ' y] ' (D,8) 



and the sum over partitions in (D.7) runs over all ordered partitions of the vector 'Y^ li N i p i 
into non-zero vectors N*/3i, a — 1, . . . , s} for s — 1, . . . , Ni, satisfying J2 a = ^ 
and 20 

EX>^X>^')> ' ( D -9) 

, a=l i=l j=l I 

for b = 1, . . . , s — 1. If Cnx^nAPi, • • • , /3 P ) given in (D.7) becomes negative then the Poincare 
polynomial vanishes. 

Based on (D.7), eq. (D.6) can be proven by collecting the terms multiplying a given 
partition on both sides of this equation and showing that they cancel. Consider for example 
a partition with s vectors of the form £)i=i Nffii for a = 1, . . . , s and let iV" = M a . Thus 
the s vectors have the form A\ + Mi/3 r , + M20 r , ... A s + M s (3 r where A\, . . . , A s are 
linear combinations of the /3j's other than (3 r . If we just focus on the contribution involving 
the U(N r ) factor then the contribution to this particular partition from the left hand side 
of (D.6) contains a factor of 

(y 2 -ir N ^H([M a ,y]\)-\ (D.10) 

a 

where we have made use of (D.7) and ignored the y~^i<j NiN ^v term since it is a common 
factor on both sides of (D.6). On the right hand side this particular partition can come from 
many different terms. In particular M a (3 r can arise as a sum of k[ a ^ copies of (3 r , fc^ copies of 
2f3 r etc. with M a = ^k^l. This in turn can arise from the Pjv 1 ,„. ) jv r _i,ixfei,ixfc 2) - ,N r+1 ,...,N p 
term if the k[ \ s satisfy fcj = k\. There is a combinatoric factor of k\\j k\ a ^\ associ- 
ated with different ways of choosing the k\ a ' copies of l/3 r out of ki copies of l/3 r . Using (D.6) 

20 Normally this relation is written as a condition on the phase of Z-^b N ap i , but we have used (3.33) 
to expres this as a condition involving the symplcctic products. 
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and (D.7) we see that the relevant factor on the right hand side is 



E E 

{fed 



{*, (a) } 



( , 2 _ ir E S E, JJ 



v-v~ 



l(yl _ y-lj 



. (D.ll) 



We can simplify this as follows. First of all using the relation N r = krf, we can express 
h in the exponents of (y 2 — 1) and y as N r — ki(l — 1) = N r — a kff\l — 1). Next 

the sum over {k{\ and the restriction Yla^i^ = h can ^ e removed if we replace all factors 
of ki by ^2 a k^ a \ Using (D.8) this allows us to express (D.ll) as 



N r „.N r 



(y 2 - iY Nr y 



nj e fe 2 -^^- i, n(^)"' 

a Aa) I \ L ' y J / 



(D.12) 



The equality of (D.10) and (D.12) is now reduced to the following identity: 

h 



[N,y}\ 



W I 



To prove (D.13), consider the equivalent identity for partition functions 

- y 2 ) 1 z l 



N 

Z ? 



E 

N=0 



[N,y}\ 



oo 



i=i 



-y 21 I 



(D.14) 



Rewriting the r.h.s. as 



I EE 

.fc=0 1=1 



[z{\ - y 2 )} 1 y 2U 
I 



I (l-z(l -y 2 )y 2 r 



(D.15) 



k=0 



we see that (D.13) follows from the known relation between the (/-deformed exponential and 
the g-deformed Pochhammer symbol with q = y 2 (see e.g. [57]), 



N 

Z 



E 



N=0 
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